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TROPICAL CURVES IN SANDPILE MODELS 


NIKITA KALININ, MIKHAIL SHKOLNIKOV 


Abstract. We consider the sandpile model on a convex part of the square lattice. We modify the maximal 
stable state by adding a grain of sand to each of n fixed points: the consequent relaxation produces pictures 
where we can see tropical curves. These curves pass through the same n fixed points and solve a version of 
the Steiner tree problem: minimization of tropical symplectic area. In order to show this, we develop several 
techniques to deal with particular integer-valued solutions of a certain Dirichlet problems. The continuous 
version of the considered relaxation reveals an interesting dynamics on polytopes. 


This is still a draft, which will be completely rewritten from the beginning. Most of the section can be 
read independently, now the paper is a set of relatively poorly connected patches. The clear statement of 
results can be found in HH. 

However we hope that the really interested reader can read this paper even in the current state. Section]^ 
contains the motivation and theorems for the rational-polygonal type of the boundary. It is outdated, the 
text is to be rewritten with a better and more concise approach. Section contains a brief introduction to 
tropical geometry. Section explains the relation with symplectic area. Section defines the dynamic on 
the polytopes. Section provides somewhat new statements about (or approach to) well-knows least action 
principle for sandpiles. Section proves the emergence of periodic patterns near the boundary. Section is 
an old version of two previous sections. 

Section 1^ shows emergence of tropical vertices, probably, the most tricky part of the paper. In Section]^ 
we prove some finiteness properties of the dynamic on polytopes. In Section]^ we start to assemble all the 
previous result in the proof of the main theorems, the case when we drop only one point. In Sections El 
|I2| we connect the emergence of tropical edges and vertices with the dynamic on polytopes, and resolve the 
degenerate cases. In Section we study the convergence questions. 

In Section 14 we define tropical series and start to consider the situation with an arbitrary boundary. In 
Section we actually deduce the situation with an arbitrary boundary from the theorems about polygonal 
boundary. Section surveys the concept of wave. 


1. Basic definitions and main statements 

1.1. Terminology. Let F be a finite graph with the set of vertices V(r) and SF be a distinguished subset 
of V(F), called the boundary of F. Vertices in V(F) \ 9F are called internal. We denote by n{v) C U(F) the 
set of vertices adjacent to v. A state is a non-negative integer function ip : V(F) —)■ No, in the following ip{v) 
stands for the number of sand grains in a vertex v. A toppling T„ at i; S V(F) \ i9F is a partially defined 
operation on the space of states. Namely, for a state (p such that (p{v) > |n(w)| — 1 we define a new state 
(p' = Ty{(p) by 

( 4>{v') — \n{v)\ if V = vP 
(p'{v') = < (p{v') -1-1 if i;' € n{v), 

otherwise. 

We say that cp' is the result of the toppling at v for a state <p. A toppling can be considered as a realization 
of privatization law: if a vertex v has at least n{v) grains of sand, then all the neighbors of v can take of v 
by one grain of sand each. Note that we prohibit to apply toppling at the vertices in the boundary of F. 
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A state (j) is called stable if (j){v) < |n(i;)| for all v £ V(r) \ dF. Thus, by definition, topplings can be 
applied only to non-stable states. It is easy to see that if each connected component of T contains a vertex 
from 9r, then after a finite number of topplings the state becomes stable. Indeed, the total amount of sand 
in I/(r) \ 9r is finite and topplings at vertices adjacent to the boundary decrease this amount. The process 
of applying topplings while it is possible is called the relaxation. This version of the abelian sandpile model 
was defined in [T]. 

The result of relaxation doesn’t depend on the order of topplings [7], that is why this model is called 
abelian. For a modern survey about sandpiles we refer the reader to |18| . Traditionally, the vertices in ST are 
called sinks but for our purposes we find more appropriate the word “boundary”. 


1.2. Explanation of results. The maximal stable state (j)o is defined by (l>o{v) = n{v) — 1 for every u € r\9r 
and (/'o('*^) = 0 for u € ST. We are interested in the result of the relaxation for states of a special type on a 
large graph which is a finite piece of the standard square lattice. Such a state differs from a maximal stable 

state by adding an extra grain of sand to vertices ui,U 2 ,... ,u„, i.e. = 4'o(.Vi) + 1 and 4>(v) = 4'o{v) 

for all other vertices. 

To illustrate the main results, we consider now a very basic example (we learn in from 0) where y(r) = 
[0, nZ^, and ST is the intersection of the boundary of the square [0,A^]^ with the lattice. Two vertices 
are connected by an edge iff the distance between them is equal to one. For such a graph the maximal stable 
state 00 is equal to 3 at each internal vertex. We define a new state 0 to be equal to 0o everywhere except 
for one point p where 0 is equal to 4. The result of the relaxation for 0 is shown on Figure We relax 0 
by doing topplings by generations, i.e. the first generation is the point p, the second generation is the set 
of vertices which became unstable after the toppling at p, the third generation is the set of vertices which 
became unstable after doing topplings at all the vertices in the second generation, etc. 

After the relaxation of this configuration we have less than three units of sand on the four segments of the 
diagonals and the sides of a small square inside, except the point p, where we still have three grains. At all 
other vertices we have three units of sand. Similar pictures appear in a general situation as explained below. 

Now we fix the notation that we use throughout this paper. Let F be a finite full subgraph of a standard 
planar lattice graph, bounded by a polygoij^ 0F be the vertices near the boundary of this polygon. Then 
we take the maximal stable state 0o and add sand grains to the points pi,... ,pn, this gives us a state 0°. 
We denote result of the relaxation of 0° by 0®"'^. On the other hand, we can drop sand to 0o into the points 
pi,... ,Pn one by one and relax after each addition. Because of the Abelian property, this gives 0®"'’* too. 

The graphic representation of 0®"'^, in general, will be very similar with the above example: vertices with 
an amount of sand less then three will be arranged along a graph (see Figure]^. Moreover, this graph has a 
natural structure of a plane tropical curve. This paper, mainly, is devoted to the proof of this result. 

In order to recognize a graph in the picture we need to assume that the number n of the points, for which 
we add extra grains, is small with respect to the size of the graph. In a stable state 0, vertices v with 
0(u) < |n-(i;)| — 1 will be called colored. We prove that the number of colored vertices in 0®”‘^ is proportional 
to the product of n and the length of the boundary of F. Since the picture appears to be scaling invariant 
and extremal in a certain sense, in the limit we have a one dimensional object, which reveals to be a tropical 
curve. This is formalized and summarized in the following subsection. 

We learned this model and pictures from [smsiis], where this problem and its properties are approached 
from the physics perpective. This paper, on the contrary, is devoted to the mathematic aspects of the 
obtained pictures. We state precise theorems and prove them via studying the limiting behaviour of the 
relaxation process. As distinct from mainstream approach started in [T], we do not study the algebraic 
properties of recurrent configurations, but consider the scaling limit of certain configurations. In some sense, 
the convergence results of our paper are analogous to the results in m and study of the limiting properties 
can be compared with [T^]. But, since the setups and techniques are quite different, this analogy is very 


limited. For the generalisations and restrictions of hypothesis, see Section 17 


1.3. Formalization. In our setup A is a lattice polygon, and pi,...,Pn are distinct points in A°, the interior 
of A. Now, let A be a positive integer. Consider a graph F^r, a full subgraph of with a set of vertices 
equal to {N • A) n Z^. Let 5F be the set of vertices of F of degree less than 4. Consider a collection of vertices 
Pi ,...,P^ € F(r 7 v) defined by taking the coordinate-wise integer part of the points N ■ pi,...,N ■ pn in A • A. 


^Later we consider an arbitrary boundary 
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Figure 1. Snapshots during the relaxation for the state (j) = 3 on a square after adding an 
extra grain at one point p (the big grey point). Black rounds represent v with 4>{v) > 4, 
black squares (which are arranged along the vertical and horizontal edges on the final picture) 
represent the value of sand equal to 2, white rounds (arranged along diagonals on the final 
picture) are 1, and whites cells are 3. Rare cells with zero grains are marked as crosses, one 
can see them during the relaxation on the vertical and horizontal lines through p. The value 
of the final state at p is 3. 


We define a non-stable state on T^ to be equal 4 at the points Pf and equal 3 at all other internal 

points. Consider the result of the relaxation process for a new stable state We define an exceptional 
(colored) set of vertices Cn in the relaxation as a set of all such vertices where is not equal to 3, i.e. 
Cn = 

Theorem 1. The sequence of sets has a limit C in A in the Hausdorff sense. Let C be the closure of 

C\dA. Then, C is a finite part of a tropical curve. Moreover, C passes through the points pi,... and the 
endpoints of C, i.e. C C dA are exactly the vertices of A. 

Definition 1. The curve C = C(A;pi,... ,p„) is called the limiting curve for a configuration of points 
pi,...,p„ G A. 

For us, the most exciting thing is that C has a natural structure of a finite part of a plane tropical curve. A 
plane tropical curve is a planar graph, whose edges e are intervals with rational slopes and prescribed positive 
integer weights w(e). The edges are allowed to be unbounded, the number of edges and vertices is finite and 
at each vertex the balancing condition is satisfied (refer to Figure]^. The latter means that for each vertex 
of a tropical curve the condition w(e)p(e) = 0 holds, where e runs over all edges incident to v and p(e) 
is the outgoing primitive vector of e. We need to recall that the outgoing primitive vector for an edge ending 
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Figure 2. Septagonal boundary. The relative position of points and the boundary remains 
the same but a scale of the grid is different. In spite of general positions of the points 
we have edges with multiplicity bigger than one. The scale of last picture is twenty times 
bigger than the scale of the first. 


at the vertex u is a shortest integer vector in the direction of the edge going out of v. For further details on 
tropical curves see Section A fundamental review of tropical curves can be found in m- 

We mentioned in the theorem that C is a finite part of a tropical curve. We call a planar graph a finite 
part of a tropical curve if it can be represented as an intersection of a planar tropical curve and some polygon 
P, such that P contains all bounded edges of this curve. Therefore the limiting curve C is a result of a 
partial cutting of all unbounded edges for a tropical curve. See an example on Figure the colored cells are 
arranged along a tropical curve; all the weights for the edges in the picture are equal to one. 

In general terms, Theorem [T] claims that near the limit the set of colored vertices will be arranged along 
a well organized finite collection of straight segments with rational slopes. Also it states that large clusters 
of non-colored vertices are convex in the limit, and for big N all non-convexity is in a small neighborhood of 
the boundary of a cluster. 

It is worth mentioning that the weights for C don’t come directly from the limit for the set of colored 
points. Nevertheless, we can extract the weights with a more delicate analysis of the relaxation process. 

Namely, instead of looking at the result of the relaxation one can count a number of topplings at each 
point. For each state on F there exist a sequence of states such that the state 0®"'^ is 
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stable and (/>* is a result of a toppling at a vertex of the state ^. As we mentioned before the state 0®"'^ 
is a result of relaxation for (jp and doesn’t depend on the sequence of intermediate steps. 

Definition 2. Let Toppl^o(y), the toppling function, be the number of topplings at v during the relaxation 
process 

The function Toppl^o depends only on and doesn’t depend on the choice of the sequence 
and the argument is the same as for the uniqueness of the relaxation. Furthermore, Toppl^o is minimal in 
certain sense, see Section [7T| 

In our setup we define such functions Toppl^^ on l/(rAr) for each N. It appears that Toppl^j,^ behaves 
well while N grows, namely, Toppl^^^ is a piecewise linear function almost everywhere. We will see that this 
linear function provides a natural choice for multiplicities for C. Finally, we will see that C enjoys not only 
the standard balancing condition at vertices in the interior of A but also certain analog of balancing at the 
boundary. This and the previous observations are made precise in the following theorem. 

Theorem 2. A sequence of functions Fn : A —)■ K given by 

FN{x,y) = ^Toppl^^{\Nx\,[Ny\) 

uniformly converges to a continuous function F'(A;pi,... ,p„) = F\ A —> K. This function is concave, piece- 
wise linear with integer slopes and vanishes at the boundary of the polygon. Then, C is the locus where AIa 
is not smooth, and thus C coincides with a finite part of a tropical curve as a set. 

As explained in the Section [7.1[ Theoremimplies TheoremIn other words, A is a tropical polynomial 
and C is its non-smoothness locus in A. Again, see Section for more details on tropical curves and 
polynomials. 

Definition 3. We call the function F(A;pi,... ,pn) the limiting toppling function for the configuration 
Pi,... ,Pn of points in A. 


Hence, a collection of points in the polygon gives rise to a curve C, and it is natural to ask what is special 
about this curve. It turns out that C{A]pi,... ,pn) minimizes the tropical symplectic area in the class of all 
curves passing through the chosen set of points {pi\ and balanced at the boundary of A. 


Definition 4 (See [3D]). The tropical symplectic area of a finite segment I with a rational slope is given by 
Area{l) = Lengthll) ■ Lengthen), where Length^—) denotes a Euclidean length and u is a primitive integer 
vector parallel to L If C" is a finite part of a tropical curve, then its symplectic area is the weighted sum of 
areas for its edges, i.e. 

Area(C') = ^ Area(e) • Weight{e). 

eeE{C') 


In Proposition it is proven that the tropical symplectic area comes as a limit for the symplectic area of 
holomorphic curves and is a deformation invariant. Usually, there is a whole smooth family of area-minimizing 
curves passing through the given collection of points, see Section 17 Therefore, minimization of area does 


not specify the curve C. Fortunately, there is another way to do that. 


Theorem 3. Let U be a set of all tropical polynomials F that vanish at dA and the curve defined by F 
passes through the points {pi}. Then the curve C{A;pi,... ,pn) has the minimal tropical symplectic area in 
the class of all curves given by polynomials in V. Additionally, F{A]pi,... ,p„) is a unique minimum for the 
functional on V given by 

F I—/ Fdxdy. 

JA 

This is known under the name Least action principle |D|. 


Definition 5. The curve given by the polynomial minimizing this integral functional will be called the 
minimal curve for the points pi,... ,Pn- 
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Thus, the theorem above states that C{A]pi,... is the minimal curve. The problem of finding such 
curve can be seen as a tropical analog for the famous Steiner tree problem: given a set of points, find the tree 
of minimal length containing all these points. 

Finally, we managed to generalize the above thoerems to the case of an arbitrary boundary. Namely, if we 
replace A by essentially any convex domain then most of the results remain true. Keeping the old notations 
it can be summarized in the following 


Theorem 4. Let C be a convex closed set with a non-empty interior such that it doesn’t contain a line 
with irrational slope. Then the sequence of sets ^Cn has a limit C in and C passes through the points 
Pi,... ,p„. A sequence of rescaled toppling functions F/v : —>• K given by 

FN{x,y) = ^Toppl^^{\Nx\,[Ny\) 


uniformly converges to a tropical series F = F: O —)■ M>o. Let C be the closure of C\dVl. Then, C is a 
tropical analytic curve defined by F. Moreover, F is the minimal tropical series with such properties. 


See Section for the discussion of tropical series and tropical analytic curves. 

The multiplicities for edges of the curve C can be found in terms of the weak sacaling limit for 
Namely, we define a sequence of functions —>■ K given by = A(3 — 0Ar([iVx])) if [IVa:] S Tat 

= 0 otherwise. 


(/'Al¬ 

and 


Theorem 5. For any smooth function ip supported on fl° 


lim 

N—¥00 


eeE 




where E is the set of all edges of C, 4 is a primitive vector of e G F and rrie is the multiplicity of e in C. 


1.4. Methods. During the first reading, we recommend to find all the definitions, remarks, propositions, 
and theorems, understand their statements - usually it is necessary to read a few lines above and below a 
statement - and only after that start reading the parts of this paper, which seem to be more interesting or 
suspicious. 

Here is the list of principal ideas and methods. 

• Looking at the relaxation process after adding a grain to (/ig we observe a sort of a ripple this point 
produces. This suggests the idea of decomposing this ripple into waves. A similar decomposition has 
also appeared in [5]. 

• Firstly, we want to know what will happen if we drop a grain to the maximal stable state. For this 
we consider what is going on near the boundary of the rational slope. Since we are interested in the 
result at the limit A —>■ oo, during this analysis we can suppose that F is just a semi-plane and we 
send waves at a point distant from the its boundary. We prove that a discrete tropical edge will detach 
from such a boundary. For that, we estimate the amount of sand, which falls to the boundary, and 
prove that the toppling function is linear near the boundary. These methods rely on the fact that a 
toppling function is the minimal solution in certain type of Dirichlet problems. 

• Secondly, we carefully studied the process how a discrete tropical edges and vertices appear from 
smoothing of a piecewise linear function. That proves that discrete tropical edges move, when we send 
waves, without changing the structure of their traceries. That is why we call such a subset of coloured 
vertices self-reproducing. 

• Thirdly, we need to globalize the above consideration. This is done via operators Gp acting on the 
space of piecewise functions. These Gp are defined in purely tropical terms, but play important role 
in the proofs about sandpiles. In fact, Gp is an “integral” phenomena while a wave is its “differential” 
version. 

• Fourthly, we need to resolve different problems of non-genericity, which appear even for generic 
situations and in the proofs. We applied different technics for that. For example, we change the 
boundary: we cut a relatively small piece of A and “freeze” it, i.e. prohibit to perform toppling 
there. Then, when the rest of the picture is stabilized, we “unfreeze” this part and prove that the 
topplings we produce now do not change the picture in the limit. Also we have to master states which 
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represents discretization of all tropical curves. Since for non-minimal curves we can not obtain their 
discretization via relaxations, we obtain them via sending waves only. 

• Fifthly, we do not need the assumption that the points pi,... ,p„ are in general position, but it is tricky 
to handle. For that, we study the dynamics of aforementioned operators Gp and prove convergence 
results about them. Then, the crucial thing is the finiteness of some product of Gp (this product 
represent the whole schema of relaxation). We prove this type of finiteness for rational pi,... ,p„ and 
then, we use the convergence results in order to localize picture near each point. Finally, we see that 
in localized picture the irrationality does not matter. 



Figure 3. Emergence of a discrete tropical edge of direction (3, 7) under the action by 
waves. It is an example of a self-reproducing pattern. 


2. A BRIEF INTRODUCTION TO PLANAR TROPICAL CURVES AND THEIR FINITE PARTS 
A tropical Laurent polynomial F in two variables is a function on which can be written as 
(1) F(x) = min{ciu + w ■ x), 

wGA 

where A is a finite subset of The numbers Cu, G K are called the coefficients of F and A is called the 
support for the coefficients of F. To emphasize the analogy with usual polynomials, sometimes we write 
P = ‘^ SujeA Cuja:’" ” using the standard tropical addition “a -I- 6” = min{a, b) and multiplication “a6” = a-\-b. 
See [4], [14] or 0 for details and motivation. 

The locus of those points where a tropical polynomial F is not smooth is a tropical curve (see M)- 




Figure 4. Examples of vertices satisfying the balancing condition. We used the balancing 
condition to define a planar tropical curve (see page|^. The notation m x {p,q) means that 
a corresponding edge has a weight m and has a primitive vector (j),q). 

A basic object that one associates with a polynomial F is its Newton polygon P{F), that is defined as a 
collection of all integer points in the convex hull of the support for coefficients of F, i.e. P{F) = I?C]Conv{A). 
Note that it is not necessary that all points in P{F) correspond to monomials in F. 

Remark 1. We always assume that, given a tropical curve G, we construct a tropical polynomial F which 
defines G such that P{F) is minimal by inclusion. 
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The polynomial F gives a subdivision for the polygon P{F) in the following way. The extended Newton 
polyhedron P{F) C is defined by 

P{F) = Conv{{w,t) G P{F) x R\t > Cw}- 

There is an obvious projection of P{F) to Taking the image under this projection of one-skeleton of P{F) 
we get a polyhedral subdivision for P{F). 

The tropical curve C defined by the polynomial F is dual to the subdivision of PiV): vertices of the curve 
correspond to faces in the subdivision, edges of the curve are in one to one correspondence with edges of the 
subdivision and faces of the curve correspond to the vertices of the subdivision. 

Note, that not all integer points of the Newton polygon are necessarily the vertices of the subdivision 
defined by F. If a point w G P{F) is not a vertex of the subdivision, i.e. it belongs to an interior of a face 
or an edge, then it means that the coefficient in the expression Q can be increased without any changes 
for the function F. In particular, this implies that a function represented by a tropical polynomial doesn’t 
determine the coefficients of the polynomial in general. 

Remark 2. So we will always assume that A = P{F) and among all the coefficients {Cu,}iugA that represent 
the same function F the one we used in Q, the definition for F, are the minimal ones. 

Note, that from now on, all the lattice points in P{F) correspond to monomials in F. 

In this paper we deal with a special type of tropical polynomials, the ones that vanish at the boundary of 
a fixed polygon A. Let T be a tropical polynomial such that F|gA = 0 and C be an intersection of tropical 
curve defined by F with A. So C is a finite part of a tropical curve. 

We want to analyze the behavior of F near the boundary. In the neighborhood of any edge e of A the 
function F can be locally written as x ^ + w{e) ■ x, where w{e) G P{F) is orthogonal to e. The integer 

vector w{e) is a multiple of a certain primitive vector, i.e. w(e) = m(e)n(e), where n(e) is a •primitive •normal 
vector to e towards the interior of A. Thus, for a finite part of a tropical curve C we constructed the function 
m on the set i?(A) of the edges of A, m = me '■ E{A) -G Z>o. 

Definition 6. The aforementioned function me is called the quasi-degree for the finite part of a tropical 
curve C C A. 

Remark 3. Note that mc(e)n(e) C P{F) for each e G F(A). We always assume that the convex hull of 
the set {me(e)n(e)} coincides with P{F), since the monomials from the outside of this convex hull don’t 
contribute to F|a. 


3. Relation to holomorphic curves 


In this section we sketch how one can translate the story about minimal tropical curves (see Definition 
to the context of classical holomorphic curves. We recommend the reader, whose interest is mostly about 
sandpiles, to skip this section during the first reading. 

Recall that an amoeba of an algebraic curve S in the algebraic torus (C*)^ is an image of S in under 
the logarithm map Log given by Log{zi, Z 2 ) = (log\zi\,log\z 2 \)- Consider a family of algebraic curves St in 
(C*)^ for t > 0. We say that the family St tropicalizes to the tropical curve C if the family of their rescaled 
amoebas LogtSt C converges to C when t tends to 00 . Here Logt simply denotes the map {logt)~^Log. 
It could seem that the topicalization C is defined only as a set. In fact, the multiplicities for the edges of C 
can be also canonically restored from the family St- 

First of all we are going to justify the name “ tropical symplectic are” that we extensively used. Suppose 
a family St tropicalizes to a tropical curve C. Let 

w = —id\og{zi) A (ilog(zi) — id\og{z 2 ) A d\og{z 2 ) 
be the symplectic form on (C*)^. 


Proposition 1. Let C be a tropicalization for St and R be a convex bounded open subset of K^. Then 


w 


47r^Hrea(C' C B) logt, 


where Bt = Log^ ^{B). 


StDBt 
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Thus the symplectic area for a tropical curve indeed can be interpreted as a main part in the asymptotic 
for symplectic areas of a family of holomorphic curves. 


Proof. The outline is the following. For a large t the rescaled amoeba Logt{St) is in a small neigborhood of 
the tropical curve C. Moreover St itself will be close to a certain lift of C to the torus (C*)^. It is performed 
by lifting each edge with a slope (p, q) to a piece of holomorphic cylinder {{z^, z’^)\z € C} translated by the 
action of the torus. This lift is called a complex tropical curve (see [13] for the details). 

Therefore, we can compute the area of St near the limit by looking at the areas of the cylinders. There 
also can be minor corrections coming from the vertices of C but the corrections are small with respect to 
logft) and so do not appear in the final statement. 

To complete the proof we need to compute the contribution from each edge in C D B. It is clear that for 
each segment in C D B the area of its lift is proportional to the length of the segment. So if we show that the 
area of the lift for the interval going from the origin to the integer vector (p, q) is equal to 47r^(p^ + 9^) logt 
then we will be done. This computation is given in the following lemma. □ 

Lemma 1. Let v = {p,q) be a primitive integer vector. Let Cf‘^ be a lift of an interval [0,u] to the torus 
(C*)^ under Logt, i.e. = {{z‘p,z‘^)\1 < \z\ < t}. Then 

/ dlog(2;i) A(ilog(zi) = —diTr^p^ logt. 

Jcr 

Proof. Let zi be rexp^icf), where r > 0 and </> S [0, 27r]. Then 

dlogzi = dlog r + icj)d(p and 


d log Zi Ad log zi = —id log r A dcjf^ 
Then the left hand side of the equality we are proving is equal to 

n 27r 

dlogr A df)^ = —diTT^p^ logt. 


□ 


Remark 4. The specific choice for uj is not crucial while it is invariant under the action of {C*Y. Indeed, 
if uj' is an arbitrary 2-form then its restriction to any holomorphic curve will not have contributions from 
pure holomorphic and anti-holomorphic parts of uj' . So we can think that uj' is a (l,l)-form. There is a 
two dimensional family of torus-invariant (1, l)-forms. Different choices for uj from this family correspond to 
coordinate dilatations on the level of tropical curves. 

Propositionsuggests us that symplectic area for tropical curves should be deformation invariant. Indeed, 
this should follow from the fact that the 2-form uj is closed. And actually, we can prove the deformation 
invariance directly. 


4 

Figure 5. In the picture we shrink a triangular cycle. Any deformation of a tropical curve 
can be decomposed into such operations or their inversions. 




Lemma 2. Consider a continuous family Cs of finite parts of tropical curves with common fixed endpoints. 
Then Area{Ca) is constant. 





10 


N. KALININ, M. SHKOLNIKOV 


Proof. Any deformation Cg locally can be decomposed into the elementary ones. An elementary deformation 
is a process of moving and shortening two edges while growing the one in the opposite direction (see Figure 


Globally this corresponds to enlarging a coefficient for a tropical polynomial. For example on Figure]^ we 
change the coefficient for the central region. 

Up to a scaling an elementary deformation simply replaces the union of segments [0,fi] an [ 0 , 112 ] by a 
single segment [0,fi + U 2 ]. Here vi and V 2 are the primitive (or appropriate multiples of primitive) vectors 
for the edges we are moving. Denote by Wi the projection of vi + V 2 on the line spanned by Vi (see Then 
after the deformation the two edges together loose 


juiljuiij + |i;2||w2| 


+ V 2 ) ■ 


Vi 

l^'ll 


+ k2|(lll +V 2 ) 


V2 

1^2 I 


kl + W2p 


of their symplectic area. On the other hand, the growing edge contributes exactly |ui +zi 2 p to the symplectic 
area of the deformed curve. □ 


W2 / 



Figure 6. Computing contributions for symplectic area. 

Lets get back to our specific case. Let the ends for a part of a tropical curve C be the vertices of a lattice 
polygon A. Suppose there is a well-defined quasidegree me for this curve. Then we can deform C to the 
union of all edges e of the polygon taken with the multiplicities mcie)- This observation together with the 
deformation-invariance proves the following lemma. 

Lemma 3. Area{C) = J2eeE{A) ^c{s)Area{e). 

We should mention that the tropical symplectic area of a tropical curve already appeared in physics under 
the name of mass of a web, where a web is a direct analog of a tropical curve. We quote [2] “On the other 
hand, we already know that shrinking an internal face of a string web corresponds to a bosonic zero mode 
(we will see that this persists in the presence of D3-branes as well in subsection 3.3), so the mass of the 
web is independent of this deformation.”, they also present this mass as a trace of some operator, and call it 
BPS-formula. “A bosonic zero mode (BZM) is an infinitesimal deformation of the web which does not change 
its mass (see [33] for a related study of the string web worldsheet).” 

Also Tension of monopolic string=Area of face and Running coupling constant=Mass of a local deformation 
mode (roughly the perimeter of the face) in BPS Spectrum of 5 Dimensional Field Theories, (p,q) Webs and 
Curve Counting 


4. Dynamics towards minimal tropical curves 

In this section we will describe certain operators Gp on the space of piece-wise linear functions, the definition 
will be done in purely tropical terms. Such an operator Gp is a limiting incarnation for the relaxation after 
dropping a grain to point p, and this relation is an important part of the proofs. It is hard to imagine that 
the remarkable properties of these operators could be ever guessed without using a sandpile interpretation. 
Moreover, with no relation to sandpiles the operators Gp are instruments in the construction of tropical curves 
with minimal tropical symplectic area. 


TROPICAL CURVES IN SANDPILE MODELS 


11 


4.1. Implicit definition of and Gp. 

Definition 7. Let i4(A,pi ,... ,pn,F) be the space of all tropical polynomials F such that 

• F{dA) = 0 

• F{p) > F{p) for any p G A 

• F is not smooth at each of the points pi. 

Lemma 4. Let A be a lattice polygon, F be a tropical polynomial vanishing on dA and pi,... be a 
collection of points in A°. Then functional F’ i-A F'dxdy defined on the space i^(A,pi,... ,p„, F) has a 
unique minimum F. Moreover, F is a pointwise minimum on V, i.e. F < F' for all F' G V. 

Fi, F 2 



L 




Figure 7. Three different relative positions for graphs of tropical polynomials Fi and F 2 
restricted to a generic line L that are not smooth at a given point. A minimum of two convex 
functions having a kink at the same point also has a kink at this point. 


Proof. The only idea of the proof is that F is a semigroup with respect to the operation of taking pointwise 
minimum. If Fi and F 2 are in V then min{Fi, F 2 ) majorettes F and vanishes at the boundary of A. To show 
that min{Fi, F 2 ) is in V it is enough to show that it is not smooth at all the points pi. Consider a generic line 
L through a point pi. Restrictions of Fi and F 2 are not smooth at pi. Then from the convexity (see Figure 
it follows that the restriction of min{Fi, F 2 ) to the line is also not smooth at pi. And thus the minimum 
itself is not smooth at all the points pi,... ,p„. 

Since taking minimum does not increase the integral we can conclude that its minimizer F exists, unique 
and is given by 

F(p) = min F'ip). 

F'ev 


□ 


Definition 8. The function Gp.^^,,^p^F is this unique minimum F in lemma above. 

The lemma has a very nice consequence. It appears that the curve defined by the minimizer F always 
minimizes its symplectic area. 

Lemma 5. For any curve C defined by F' S F 

Area{C') > Area(G), 

where G is given by the polynomial F. 

Proof. Point-wise minimization implies the point-wise minimization for quasidegree. To see this, one needs 
to look at the neighbourhood of each edge of A. Now the statement follows from the deformation invariance 
for symplectic area in the form of Lemma □ 

Now, as an initial data we have a point in the interior of the lattice polygon A. We transform a tropical 
polynomial F vanishing at the boundary of A to a tropical polynomial F' also vanishing at the boundary of 
the polygon. Also, the curve G' given by F' will pass through p, i.e. the function F' will be not smooth at 
p. We will show that for any F there exist a canonical way to find such a tropical polynomial and we denote 
it by F' =: GpF. 

Definition 9. Let F(A,p, F) be the space of all tropical polynomials F such that 
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• F{dA) = 0 

• F{p') > F{p') for any p' S A 

• F is not smooth at p. 

Proposition 2. For any point p G A° and any tropical polynomial F vanishing on dA the functional 
F i-G J^F defined on the space V{A,p, F) has a unique pointwise minimum F'. 

Proof. We proceed analogously to the proof of Lemma □ 

Definition 10. We denote by GpF this unique minimum F'. 

Essentially, Gp can be seen as a result of squeezing the face $ where the point p sits, until the curve and 
p meet (see Figure [^. For more details see Section 

4.2. Recursive definition of Gp^^,,^p^. Now we start with a sequence of points pi,..,p„ and the tropical 
polynomial F vanishing at the boundary of A. We construct a certain polynomial F' = such that 

the curve it defines passes through all the given points. We will define the new function recursively. If F is 
not smooth at all the points pi, ..,pn then we take F' to be F. If it is not the case, i.e. there exist p in the 
set {pi} of the given points such that F is smooth at p, we define Hp^^,,^p^F as Hp^^,,^p^GpF. It is completely 
unclear from the very beginning that this recursive definition always terminates and the result doesn’t depend 
on the choice of the point p G {pi, .■,Pn} at each step. 

Lemma 6. Let F be a restriction of a tropical polynomial to the polygon A and pi,... ,pn be a collection 
of points in A°. Then 

Gpi,...,p^F{x) > Gp^^ .. .Gp^^F{x), 

where for every point x G A and ki = I,..., n. 


4.4 


Proof. For m = 0 the inequality is trivial. Suppose we know that the inequality is true for m — 1, i.e 
Gpj^^...^p„F > Gp^. ...Gpi^^F. It is clear that Gp^^,,,^p^F is not smooth at Pm and so using the implicit 
definition (Def. of Gp^ we finish the proof. □ 

Therefore, passing to the limit, we have > F[p^^,,,^p^F{x). From the other hand, G 

V{A,pi,... ,pn,F) and it is not smooth at all pi. Therefore, by definition Gp^^,,,^p^F{x) < F[p^^,,,^p^F{x) and 
we proved the following lemma. 


Lemma 7. The operator is well defined, does not depend on the order of points and coincides with 


pi....,p„; but, depending on the situation, we will treat Gp^^ 


We will never use the notation iJj 
the minimum of the functional, or as the limit of some product of Gp. 
enough to apply only finite number of Gp. that is proved in Theorem 


. In fact, in order to get 

Also, it appears that Gp^, 
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either as 
it is 
is a 


piece-wise linear analog for the relaxation process, and this relation is formalized in the following paragraphs. 


4.3. Back to discrete tropical curves and their limits. Now we return to the problem of describing 
the limiting curve G obtained as the limit for relaxations of 0^. Suppose that we already know the curve 
C{A;pi, ..,pn). Consider another point Pn+i in the interior of A. It is natural to ask how the new curve 
G(A;pi, ..,pn+i) can be derived from the old one. 

There are two very different cases. The limiting curve either passes through Pn+i or doesn’t pass through 
it. When the new point is on the curve we claim that a new limiting curve and a new limiting toppling 
function are the same as the old ones. Recall that F is the limiting toppling function. 

Lemma 8. If the point Pn+i G A° is on the curve G(A;pi,... ,pn) then 

F(A;pi, ..,pn+i) = F{A;pi, ..,p„). 

In particular, G(A;pi, ..,p„+i) = G(A;pi, ..,p„). 

Proof. □ 


In the case when Pn+i is not on G the situation is more complicated. It appears that the operator Gpj^.._p,^^^ 
is precisely the one transforming the old curve to the new one. 
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Proposition 3. For any collection of distinct points Pi, ■ ■ ■ ,Pn+i G 

■ ■ • jPn+l) — J i^(A, Pi, . . . ,Pn) 

Proof. See Proposition]^ □ 

Propositiontogether with Lemmas [7| and the construction of give the following remark. 

Remark 5. F(A;pi,... ,p„) = Gp^^,._p^ (0), where 0 denotes the function on A which gives zero identically. 
Theorem 6. The curve given by the function Gp^^,,^p^{0) is minimal for the collection of points pi, .■,Pn- 
Proof. Remark and Theorem show that the operator Gp^^,,^p^ provides an entirely geometric way to 


construct a minimal curve, for details see Section 13 


□ 


4.4. Explicit description of Gp. In the remainder of the section we are going to give an intrinsic description 
for the operator Gp. This can be seen, perhaps, as a constructive version of the second part of Proposition 

m 





Figure 8. Action of Gp by shrinking the face $ where p belongs to. Firstly, t = 0, then 
t < tf), and finally t = to. Note that combinatorics of the new curve can change. 


Let us take a point p in the interior of A. The construction of the new polynomial Gp will depend on 
position of p with respect to G. 

The definition of Gp implies that if the curve G passes through a point p then GpF = F. 

Suppose that p is in the complement to G in A. Then we have two more cases. Let $ be the connected 
component, i.e. the face of A\G to which p belongs. The face can be adjacent to the boundary of A or not 
adjacent to it. 

If 4) doesn’t touch the edges of A then the procedure is quite straightforward. Consider a vector w S P{F) 
that corresponds to 4>. If we are increasing the coefficient Cyj we squeeze the face 4) (see Figure This 
produces the family of polynomials Ft such that Ft differs from F only by a coefhcient before the monomial 
“a;’"” which is defined as c^, + t. We take Gp = Ft^ such that to > 0 is the smallest number for which the 
curve given by Ftg passes through p. 

In the case when 4) has common edges with the polygon A we cannot do exactly the same thing because 
we will violate the vanishing at the boundary. To avoid this we increase the quasi-degree for the curve G and, 
if it is required , we extend the Newton polygon of F. 

In fact, if EIa is not identically zero and G H A° is not empty, then the set 4) C 9A consists of at most 
one edge. So we suppose for now that G is not empty and 4> H dA = e. We can define a new curve G', which 
coincides as a set with G incide A, and a new polynomial F' by modifying the old ones, i.e G' = G U e and 

F'{x) = min{F{x), {mc{e) + l)(n(e) • x + Og)), 

where Og is chosen in such a way that the function x i—>■ n(e) ■ x + Oe vanishes on the edge e. Note that F' 
coincides with F inside A. So we have essentially reduced the situation to the previous case. We can start 
increasing the coefficient for w = mc(e)n(e) and so forth, until a curve crosses the point p, thereby defining 
GpF := GpF'. 

The only thing that really matters in this approach is the accurate choice of the polynomial representation 
for F' (for F in the previous case). As we explained in Section]^ we should assume that the coefficients for F' 
which are not the vertices of the subdivision for the Newton polygon P{F') are chosen to be smallest possible. 
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Afterwards we will give another constructive description for Gp free of these technicalities (see Proposition 

181 ). 



Figure 9. First row shows how minimal curves for one point on a pentagon boundary 
depend on the positions of the point. The second raw shows a dual decomposition for their 
Newton polygon. 

The last case is very special, it is the case when <1)° = A°. The curve C itself is not necessarily empty but 
it can be the a union of some edges in the boundary of A. Following the analogy with the previous case, we 
define a new curve by completing C to the whole boundary, i.e. C = 9A, and in the same way we construct 
a new function F' by taking the minimum of all linear functions n(e) ■ x + for e € F(A). Then GpF is 
equal to GpF'. And as before GpF' is derived from F' by increasing the coefhcient cq. The dependence of the 
result on the position of point p is shown in Figure 

Lemma 9. The above paragraphs give direct construction for Gp from Definition 

Proof. Indeed, GpF is the point-wise minimum and during this construction we increase only one coefficient 
of F as less as possible to have p in the non-smooth locus. □ 

5. Least action principle for locally finite relaxations 

In this section we recall some basic facts on abelian sandpiles. Sandpiles can be defined on arbitrary 
(multi-)graphs, see for the details. Most of the following statements and concepts are well-known in slightly 
different contexts, especially for abstract finite graphs. 

Consider a set D C such that fl is closed and has a non-empty interior. We define a graph F = r(D) 
with the set of vertices equal to D n Two vertices vi and V 2 are connected by an edge if they are adjacent 
in the lattice Z^, i.e. if |ui — U 2 I = 1. So, the edges of T are not oriented and T is a graph without loops and 
multi-edges. We allow D to be non-compact, therefore T can be an infinite subgraph of Z^. 

We will assume that T is non-empty and connected graph. In what follows, u G T means that u is a vertex 
of r. When we say about a function on T, it is always a function on the vertices of T. For each u S F we 
denote by n{v) the set of all adjacent vertices (neighbors) of u in Z^, so |n(u)| = 4 is the degree of a vertex v. 

Definition 11. The boundary 9r of the graph T is the set of all the vertices of Z^ \ T adjacent to T i.e. 
9r = {u S Z^ \ r : n{v) C T 7 ^ 0 }. 

5.1. Topplings. A state <() is a non-negative integer-valued function on the set of the vertices of T. States 
represent distributions of sand on the graph T. 

For each vertex u S F we define the following operator Ty on states called a toppling at v: 

{ (fiw) — 4 if u = w, 

(l){w) -I- 1 if u e n{w), 

4i{w) otherwise, 

where ^ is a state and w G F. Note that we never do topplings outside F. Following the terminology of , 
for a given state (f a toppling Ty is called legal if Tytj) is again a state. In other words, the toppling Ty is legal 


30 
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for (p if Ty(j)(v) > 0 or, equivalently, (piv) > 4. The vertex v is called unstable (resp. stable) for p, if p{v) > 4 
(resp. p{v) < 4). 

The passage from p to Typ can be interpreted as a redistribution of four grains of sand at v to its neighbors. 
By definition of the toppling operator, during this redistribution the amount of lost sand is equal to |n(p)n(9r|. 

5.2. Relaxations. Consider a state p. If there are no legal topplings for p, i.e. any v e F is stable for p, 
then p is called a stable state. 

If p is not stable, then there exist a vertex vi such that the toppling Ty^ is legal for p. Applying Ty^ we 
obtain a new state pi = Ty^p. If pi is not stable, i.e. pi{v 2 ) > 4 for some V 2 G F, and we can again obtain a 
new state p 2 = Ty^pi. Repeating this procedure we obtain a sequence of states {pi\. This process is called a 
relaxation for the state p. In general, there are many ways to perform a relaxation for a given state p. 

Suppose that a relaxation terminates after a finite number of steps (which is always the case when O is 
a compact set). Then it is well known that the resulting stable state depends only on p and doesn’t 
depend on the order of legal topplings performed during the relaxation. The stable state is called the 
result of the relaxation. 

In the case of infinite F, a sequence of legal topplings can be infinite. In fact, for two relaxations the resulting 
states (if they are well defined) can be different. Thus, the questions what is the result of a relaxation of a 
given state and whether it is well defined are more delicate in this case. 

Definition 12. A relaxation {pi] 

p = po ^ Pi ^ ^ pi ^ . 

where pi = Ty.pi-i and € F is called locally-finite if for each w € F the sequence pi{v) stabilizes. 

In other words, for each v € F there is only a finite number of topplings performed at v during the relaxation 
{Pi}. Counting the numbers of topplings at a vertex of F we define a toppling function (odometer) H given 
by H{v) = |{z : Vi = i;}|. We extend H to 9F by zero. This is natural since we never do topplings at the 
vertices of i9F. 

For a locally-finite relaxation sequence we define a resulting state 

Pooiv) = lim pi{v) = po + AiJ 
2—>-00 

that is determined by the toppling function H. Here A: —>■ if denotes the usual discrete Laplace 

operator 

Y, {F{w) - F{v)) = -AFiv) + E 

w^n{v) w^n{v) 

Definition 13. A locally-finite relaxation sequence is called stabilizing if the resulting state poo is stable. 

We give a simple criteria for the existence of a stabilizing relaxation sequence and uniqueness of the 
resulting toppling function. 

Proposition 4. For a state p there exist a stabilizing relaxation if and only if there exist a function F : F —> 
Z>o such that p + AF < 4. Furthermore, in this case the toppling function for the stabilizing relaxation is 
pointwise maximal among toppling functions for all relaxations. 

A state p for which there exist a stabilizing relaxation will be called relaxable. The resulting state p + AH 
is the result of the relaxation for p. 

Example 1. Note that not every state is relaxable. For example, = M.‘^,p{v) = 4 for each v G Z^ is not 
relaxable. 

Indeed, suppose that there exists a locally finite relaxation, we look at its toppling function H. Let v be a 
point where H attains its minimal value. Counting of numbers of incoming and outgoing grains at v shows 
that p^^‘^{v) = 4 — 4:H{v) + ^yj^y H{w) > 4 , which is impossible because must be at most 3 at every 
point. 

Proposition 5. Note, that if Z^ is not contained in O and the set {v|(()(u) > 4} is finite, then p is relaxable. 

We prove this proposition later in order to demonstrate several technics we develop (duality of sand for 
for n = Z^ \ {(0, 0)} and locally finite relaxations merge). 
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Lemma 10. If there exist a function F: F —>■ Z>o such that (j) + AF < 4, then for each relaxation, its 
toppling function H is pointwise less than F 


Proof. Assume that such F exists. Consider a locally-finite relaxation 


</* — </’0 (/>i , 


where (pi = Ty.pi_i, and its toppling function H. We show that H < F. 

Denote by a sequence of intermediate toppling functions, i.e. 

(2) Hm{v) = \{i ■. i <m and Vi = u}|. 

Clearly, F^+i > Hm, (pm = P + AF^ and 

H{v) = Hao{v)= lim Hmiv). 

m—^oo 

The function F is non-negative so F > Fq. Assume that F > F^- To prove the inequality F > F^+i it is 
enough to show that F(ti„+i) > Hm+li^ra+l) = Hra{Vm+l) + 1- 

We have pmi'Vm+i) > 4 because is a legal toppling for (pm- Rewriting this and applying the inductive 

assumption we get 

4 ^ (pmiy-m+l) — Pig^m+l) T — 

— 4F^(Um+i) “b ^ ^ Fi 

< p{Vm+l) - ^Hm(Vm+l) + ^ F{w) = 

< 4 -b 4F(um+i) - 4F^(w^+i). 

Therefore, F(um-i-i) — Hm{vm+i) > 0. Since F^ and F are integer valued 

n Vm+l) ^ dUmiym+l] “b 1- 


□ 


Proof of Proposition^ Now we will construct a stabilizing relaxation for p. Consider a set R of all locally- 
finite relaxations of a state p. Let F*" be the toppling functions for a relaxation r G R. Note that by definition 
every finite sequence of legal and gives a locally-finite relaxation sequence. Thus the set R is non-empty, in 
particular 0 G R. Consider a function F G Z>g given by 

R{v) = maxF’’(u). 
reR 

This function is well-dehned and F < F. We claim that H is a toppling function for a stabilizing relaxation. 
Lemma 11. For each two finite relaxations ri,r 2 G R the function 

H(v) = max(F’'^ (u), F*"^ (u)) 
is a toppling function of some relaxation r G R. 

Proof. Consider intermediate toppling functions Flf ^, F[^ as in ([^ . We construct a sequence of functions 
Fi : r —>■ Z>o. Let Ni,N 2 be the numbers of steps in the relaxations ri,r 2 respectively. Let Hi = for 
i = 1,2 ... ,Ni and Hi — max(FjVi, Hf^) for i = Ai -b 1, Ai -b 2,..., Ai -b N 2 . It is left to the reader to verify 
that {p + AHi){v) > 0 for u S r and i = 1,..., Ni + N 2 . Therefore, if we erase from the sequence Fi, F 2 ,... 
repetitions (it is possible that Hi = F^+i on F, in this case we erase F^+i), the resulting sequence represents 
a relaxation process. We call the resulting relaxation r = m{ri,r 2 ) the merge of ri,r 2 . Note that m(ri,r 2 ) 
coincides with ri at the beginning by construction. □ 

Remark 6. The same proof works for two non necessary finite locally finite relaxations ri,r 2 . Indeed, 
consider a sequence of compact sets Di C fl 2 C ..., IJ = D. Then, consider Ni G N such that after 

Ni steps ri,r 2 has no topplings at vertices in F n Let be the merge of ri,r 2 up to step Ni. Then, 
do the same for all Di, i.e. is the merge of F and r'i,r 2 up to step Ni,i = 1,2,.... Note, that 

... coincide on F H therefore the limit r°° = limi_>.oo F is well-defined locally-hnite relaxation, 

and H’^°° (v) = max(F’’i (u), F*"^ (u)) for all u S F. 
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The set of all the vertices in T is at most countable so we can represent it as a sequence of vertices 
r = {wi,... ,Wk, ■ ■ ■}■ For each Wk we find a finite relaxation with ik steps such that Hl^{wk) = 
T-L(wk)- We merge these relaxations and obtain a relaxation corresponding to %. Namely, the sequence 
ri, TO(ri, r 2 ), m(m(ri, r 2 ), ra),... is a sequence of tails of a relaxations whose toppling function is %. Note 
that (f) + AH is a stable state, since otherwise one could increase the value of H at a certain point. 

The other implication is trivial. If there exist a stabilizing relaxation then its toppling function is integer¬ 
valued non-negative and vanishes at the boundary. Since the resulting state is stable we have (p + AH < 4. 

Finally, if we have two toppling functions Hi and H 2 of certain stabilizing relaxations than (p -|- AHi < 4, 
thus Hi > H 2 and H 2 > Hi, so Hi = H 2 . □ 

Corollary 1. Any relaxation sequence of a relaxable state is locally-Hnite. 

This was not evident from the beginning. A priori, one could imagine that there is a locally finite relaxation 
which gives a stable state, but there are other, non-locally finite, relaxations. 

5.3. Waves. In this section we introduce a classical concept of waves in sandpiles. 

Definition 14. The basic operator S'„ of adding one grain of sand at some vertex v and then relaxing . 

We will later show that Sy essentially decomposes into a sequence of waves. 

Consider a stable state (p and a vertex i; e F. Let (po be the result of adding an extra grain of sand to a 
vertex v, i.e. po = (p + Sy. The result of the relaxation process for pQ can differ significantly, for different v, 
we classify the possibilities. 

• The state po could be stable. This is the case when p{v) < 3 and the result of the relaxation is simply 

Po¬ 
rn Then, po could become stable after a single toppling. This would mean that p{v) = 3 and Syp = Typo- 
In particular, p{w) < 3 for any w G n{v). 

• Finally, the relaxation could go in a nontrivial way and require more than one toppling, i.e. an 
avalanche of topplings. 

The last case appears when v belongs to a cluster of vertices at which p takes 3 as its value. We call such 
clusters territories. Here is a formal definition of this concept. 

Definition 15. A subset T of F is called a territory of a stable state p if 

• T is connected as a subset of the graph F 

• p{v) = 3 for oXIvGT 

• p\w) < 3 for all w G F\T such that n{w) n T ^ 0 

• T has at least two elements 

The first three conditions imply that T is a connected component of the subgraph where p = 3. The last 
condition is necessary to avoid trivialities (second case in the above list, when Syp = Typo) and to give a 
correct definition of waves. 

Suppose that a vertex v belongs to a territory T of a stable state p. For now we also suppose that po = p+Sy 
is relaxable. The relaxation for the state po requires more than one toppling. Since the state p is stable, Ty 
will be the first toppling in any stabilizing relaxation of po- Denote by pi the result of applying T„ to pQ. 
After we apply a toppling at each non-stable vertex from n{v) PlT we obtain a state p 2 - Note that p 2 {v) > 0 
and p 2 is relaxable since it is a result of a partial relaxation for po- Thus we can define a state po = P 2 ~ '^v 
and by Proposition (()3 is relaxable because p2 = po + 5 y is relaxable. 

Definition 16. Let Wyp be the result of the relaxation for po- We call Wyp the result of sending a wave 
from V on p. 

Remark 7. Note that Syp = SyWyp if ?; is a vertex in some territory. Moreover, since each wave sent from 
V produces at least one toppling at v and the state p + 5y is relaxable by the assumption one can send only 
a finite number m of waves until v is not in the territory of Wt^p. Thus, Syp can be written as SyW^p and 
SyW^p differs from W)J^p by adding a grain at v and performing at most one toppling at v. In this approach, 
Sy represents an avalanche of topplings which is splitter into a sequence of “infinitesimal” waves Wy. 

Lemma 12. If vi and V 2 are in the same territory of p then Wy^^p = Wy^p. 
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Proof. It is enough to check this statement for the case if V 2 S n{vi), which is a straightforward calculation. □ 

This allows us to define W-j-f), the result of sending a wave from a territory T, as if i; is a vertex in 
a territory T of a state 4>- 

Proposition 6. Let v, v' be two vertices in a connected component of T, (f is any state. Then, there exists a 
number N such that there exists a relaxation of (f + N5y, such that we perform a toppling at v' during this 
relaxation. 

Proof. Consider a path I in T from v to v'. Let L be the length of 1. Note, that if we add 4^+^ grains to u, 
then we can perform a relaxation such that we have at least topplings at v, at least 4^“^ topplings at the 
neighbor of v in /, at least 4^“^ topplings at the next neighbor in 1. etc, at least 1 toppling at v'. □ 

Proof of Proposition^ It is enough to consider the case when H contains all the lattice points except (0,0), 
r = \ {(0,0)}. It follows from Remark § that it is enough to consider (p, such that ^ = 3 at all vertices 

except one vertex v, where (p{v) = 4, and prove that (p is relaxable. So, we need to prove that there exists N 
such that (j) has at most 2 grains at v. Consider Wffp for some n S N. During this process, some number 
A of grains have fallen to (0, 0). We add (0,0) to T, consider Wffcj) and add Al + 3 grains to (0,0), this gives 
a state ip. Note, that the relaxation of ip gives the state where we have 3 grains of sand everywhere. Now, by 
Proposition there exists such big N that if A > IV then during the relaxation of ip we add grains to v. □ 

Proposition 7. If T = Z^ \ {(0,0)} and the set {u|^(u) > 4} is infinite, then (p is not relaxable. 


6. Pattern formation on a half-plane with rational slope 


Consider the lattice graph Z^ with a trivial boundary i9Z^ = 0. Consider a line La with a slope a G PQ^ 
such that La passes through the origin. Let (p,g) be a primitive vector in the direction of La, i.e p and q 
coprime integers and [q : p] = a. Denote hy F = Fp^q a linear function T’: > K given by F(x, y) = qx — py. 

The zero locus of F is La. Since p and q are coprime we can choose two coprime integers (p',q') such that 
F{p',q') = l. 


6.1. Formulation of the main results. Here we prove that there exists a unique periodic soliton state 
with any given rational slope. 

Definition 17. A sand soliton in the direction [p : q] = a G PQ^ is a stable state p on iP such that 

• p{{vi,V2)) = Pii-vi, -V2)) and p{{vi,V2)) = </>((ui +p,V2+ q)), for any {vi,V2) G I? 

• the set {u G r|(/>(ii) < 3} is contained in a finite tubular neighborhood of La- Note, that the set 
{u G r|())(i') = 3} contains exactly two infinite connected components which we denote by F and —T, 
such that the restriction of F to T (resp. —T) is positive (resp. negative) 

• W±rp{{vi,V 2 )) = 0((ti=Fp', V 2 Tq')), for any (ui, T 2 ) € Z^, i.e. a wave acts on the set {u G r|<)>(u) < 3} 
by the translation by the vector {p,q). 

Theorem 7. For each a G PQ^ there exist a sand soliton pa of the direction a. 


Proposition 8. All the solitons in a direction a differ by a translation. 

The pattern from Theorem naturally appears during certain relaxations on convex domains. The most 
basic example of these relaxation is the following. 

Consider P = r(D) a graph with boundary (see Definition 11), where 

D = {(x,y) : F(x,y) > 0}. 

The set D is a closed half-plane such that 9D = La. Let v be an interior vertex of T and p be the state on T 
given hy p= (3) -I- 5^. 

Lemma 13. Consider two functions Fi and F 2 defined on a subset of 1? , such that AFi < 0, AF 2 < 0. Then 

Amin(Fi,F2) < 0. 

Proof. Consider a vertex u G T. Let Fi(v) < F 2 {v). Then, Amin(Fi, F 2 )(u) < AFi(u) <0. □ 


Lemma 14. The state p is relaxable and the toppling function Hy of the relaxation is majoreted by Fy{w) = 
Taiii{F{w),F{v)),w G T. 
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Proof. Indeed, AF„ < 0 by Lemma 13 and AFy{v) < 0 by the construction. Thus, the statement follows 
from Proposition 1^ since 4> + AF„ < 4 and is positive integer-valued function on P. □ 

Lemma [l4| implies that one can speak oi cj) = Sy{3), that is the relaxation for the state (p, where (3) denotes 
the maximal stable state. The following theorem describes an approximate shape for this state in the case 
when V is sufhciently far from the boundary. In particular we observe certain stabilization phenomenon and 
emergents of a soliton pa- 

Theorem 8. If the distance from v to is greater then da than the defect set {w G r|(/)(w) < 3} is contained 
in the union of two tubular neighborhoods i?i and B 2 of La and (La + v) of radius Tq,, where Va and da are 
positive absolute constants depending on a. The restriction of </> = S'„3 to Bi doesn’t depend on v and is 
invariant with respect to the translation by (p,q). Moreover, the restrictions oi p — 6y to B 2 and pa to Bi 
coincide up to a translation. 

6.2. Proofs. We start with Theorem]^ and then using it we deduce Theorem 

Proof of Theorem^ By Remarkit is required to send a finite nu mbe r m of waves from v before v is outside 


a territory T of IT™(3), where p = (3) + Sy. Note that by Lemma 12 sending a wave from v on each W^(3), 
for fc < m, is the same as sending a wave from any point with the same or a greater distance from La- 
Therefore, each state W^(3), for k < m, is invariant with respect to the translation by {p,q). Here, by a 
abuse of notation we denoted by T the territories for different states W^{3). In fact, T is destinguished from 
all other territories by the condition that the restriction of F to T is unbounded. We interpret the operator 
Wj- as sending a wave from infinity. 

We want to estimate the amount of lost sand after sending a number of waves. The naive count suggests us 
that the number of lost sand should be infinite since there are infinetly many sinks and the state is periodic. 
To avoid these problem we can take a quotient of the whole picture by (Np, Nq) , where TV S N can be chosen 
arbitrarely. Let T' and 9r' will be the result of taking this quotient, Wy3', F' and denote the push-forward 
of W^{3)^ F and F„ (see Lemma 14) to T'. One can think of T' as a half-infinite cylinder. 


Note that the number of topplings at the verteces adjacent to dT is bounded by F„ which is itself bounded 
by F (see the beginig of this section). This means that after sending a sufficient number of waves M, the 
waves stop reaching the boundary. More formally, if we denote by Hk a toppling function corresponding to 
the partial relaxtion for (3) -I- Sy after sending k waves from v, for any ru S T an sequence integer F[k(w) is 
non-decreasing and bounded by F(w), and so Hk{w) stabilizes. In particular the state W!f3 stabilizes at any 
vertex. Denote by a state which is obtained as a result of a pointwise stabilization for WJf{3), i.e. 

Wf’(3)(w) = lim Wif{3)(w). 

A;—>-oo 

The uniform bound by F would also imply that cardinality of a deffect set {Wyp' < 3} is bounded if the 
graph r' was finite which is not the case. But we can use a refined version of the same argument. 

Lemma 15. Let be a push-forward of Hk to Lb Then 


AF' = 


Hi 


I dr' 


Note that the first integral — AH(, is well defined since is constant sufficiently far from Lo 
integral counts defect verteces with positive multiplicities. 


Proof. 

The right hand side in the equality of Lemma 
number of defect vertices. 


This 


□ 


15 


is bounded by F' for any fc. Thus, has only finite 


The definition and the structure of Wf?3 allows us to give the first estimates on the constants Tq 


and da- 

Let r be the maximum of distances from La to all defect vertices of Wf?3. We will take Tq, to be greater 
than r. We have already shown that there exist such d that for all v the states Sy3 coincide with in a 
neighbourhood of La of radius r. We will take da to be at least d. 

Now we can destinguish two types of defect verteces on W^3 for fc large enough. The ones that are stabilized 
and those that “move”. Moreover, we can destinguish a new territory 71-, which separates satibilized and 
moving defect verteces (see ). 
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We are going to use the following fundamental theorem (see 0 ). 

Theorem 9. Let i? be a positive real number. Consider a real-valued function / dehned on a set {v € : 

|u| < R}. Suppose that / is descrete harmonic and non-negative. Then 

1/(0,0) - /(1,0)1 < ^ maxf(v). 

R |D|<i? 


Morally, this theorem can be seen as an estimate on a derivative of a harmonic function. We call dxf{x, y) = 
f{x -f 1, y) — /(cc, y) a descrite derivative in the ai-direction. A derivative in the direction of y is dehned in 
the similar way. 

We want to apply this theorem to the restriction of Hk to the territory 7/. Consider a point w G r° such 
that dist(w, La) > r + r' for some r' > 0. It is clear that for k large enough w is arbitrary far from the 
set of moving deffect verteces of W^3. In particular w is contained in Tk together with a disc of raduius 
r'. Applying Theorem]^ we have d,Hk{w) < Hk{w)/r', where d, means either dx or dy. Using the theorem 
second time and applying trivial estimates we get 


1 f) n ( \ ^ Hk{w) Fjw) _ dist{w,La)^p^ + q^ {r + r')^p^ + q^ 

* • k 

Clearly, there exist rg > 0 such that for any ri > rg 


(r -I- 

ri(r-i - 1) 


Therefore, if r' > rg then di,d,Hk{w) = 0. Thus, the restriction of Hk to Tk is linear. 


□ 


7. The criticality of the toppling function and the structure of self-reproducing 

PATTERNS 


7.1. The criticality of the toppling function. Given a state we denote by hy the number of topplings 
at the vertex v during a relaxation process of so hy = Toppl^oiy) in the terminology of the first chapter. 
The inequality 

(3) 0 < (l3°{v) - \n{v)\■ hy + ^ < |n(u)| - 1 = 0o(^’) 

w^n{v) 

follows from counting the number of incoming and outgoing grains at v. Recall that (j)g is the maximal stable 
state, i.e. (I)g{v) = |n(u)| — 1 for all vertices v. 


Definition 18. Let Ahy = J2wen{v) — hy be the Laplacian of h. 

If = (pgiv) = \n{v)\ — 1 for a vertex v then the inequality ([^ implies 

-1 < Ahy < 0, 

and it means that hy is a superharmonic function. If (f>^{v) > |u(u)|, i.e. we added grains there, then 


(4) 


Mv) - <P°{v) _ 1 < AL < Mv)-fiv) < 0 


|n(z;)| 


|n(u)| 


Remark 8. Note that Ahy = 0 means that (j)'^{v) = This implies that the main object of our study, 

i.e. the set of vertices where differs from (jig, has other description: 

{u|(()o(^’) T = {v\Ahy < 0}. 

In the limit, the function h becomes piecewise linear with integer slopes (Theorem]^. Hence the limit of 
the set where Ah < 0 is exactly a non-smooth set of such a function, i.e. this limit is a tropical curve by 
definition. Note that we have Ahy < 0 at a vertex v where we have added sand (i.e. cj)°{v) = (j)g{v) -I- 1), 
because of Q, so u belongs to the tropical curve. Then, Theoremimplies Theorem]^ 

We fix an order of topplings during the relaxation of (/)°, i.e. a sequence • • • —of states. 

We consider functions /i^, ..., where hi = F{v) is the number of topplings at v among first i topplings 

during the relaxation. Note that the possibility of applying the toppling to ((>* at v is equivalent to A/i(, > 0. 
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We state a new problem. Find a function H : T ^ Z>o,i?|ar = Oj which has AH < 0 at the internal 
vertices, and AH < 0 at the fixed different vertices pi,P 2 t ■ ■ ,Pn- This initial data correspond to the initial 
state (fP which differs from (pQ by added sand grains at pi,p 2 ,... ,Pn- The sand relaxation of (j)° provides us 
some solution for the problem of finding such H. 

Proposition 9. The sand solution (the function is the pointwise minimal solution of the 

above problem. 

Proof. This follows from the fact that each of h* is smaller or equal at every point than any solution H. 
Indeed, for = 0 it is clear, < H. Let h® be less then H at every point. We will show that this implies 
the same inequality for Let differ from h® by the toppling at w. If w does not coincode with one 

of Pi then Ah\^ > 0, and the latter implies that 

\n{w)\ ■ + 1 < (-1 + ^ hl) + l< ^ Hy <\n{w)\- H.^,- 

vGn{w) vGn{w) 

Since h\y, Hy, € Z, we proved that h^y + 1 < Hy,. Similar argument works it w = pi for some i. □ 

This is a lind of minimal action minimal action principle, “sandpiles are lazy” in the terminology of Levine 
etc. 

7.2. Waves. During the relaxation process we see a sort of ripple effect, see Figure So we decompose the 
relaxation as a sequence of waves ([mils]). Composed in a different order, the waves constitute the operators 
Gp, that explains why the relaxation in the limit is a sequence of the operators Gp. 

Definition 19. Given a stable state (p, consider the full subgraph that contains all the vertices v such that 
(p{v) = \n{v)\ — 1. Consider the connected components of this subgraph. Those who contain at least two 
vertices are called territories. 

Definition 20. (ini) Let u be a vertex in some territory in a stable state. The wave at v is the following 
procedure. We add a grain to v and topple v if possible. Then, we “freeze” u, i.e. prohibit to topple it and 
perform a relaxation. After stabilisation we remove a grain from v. 

Remark 9. It is easy to see m) that after “freezing” of v and subsequent relaxation, each vertex can 
topple no more than once. Basically, a wave increase the function Toppl by one in this territory. 

Lemma 16. The adding of a grain to v and the relaxation makes the same effect as making waves at v while 
it change something (i.e. until we have (p{v) < |n(u)| — 1), then adding a grain to v. 

Proof. Left to the reader. □ 

On the video of the relaxation process (Figure]^ we will see that a wave moves the set of colored vertices 
towards v. We will prove that the set of colored vertices is made of self-reproducing traceries, and under this 
action by waves this tracery finally crawls over u, and the relaxations terminates. In the limit this crawling 
is easily recognized as the operators Gp. 

Proposition 10. The results of waves made at two different vertices in one territory coincide. 

Proof. Consider a state (p and two adjacent vertices v and w in the same territory. It is enough to prove that 
the results of waves at v and w coincide, that is straightforward. □ 

7.3. Discrete Harmonic functions. Let By{p) C 1? be the set of points in 1? with the distance less than 
r from p. 

Lemma 17. ([5], Theorem 5) Let u{p) be discrete harmonic and non-negative function in Bn(j>Q),R > 1. 
Let Pi be a neighbor of po: then |it(pi) — u(po)| < Gu{po)/R, where G is an absolute constant. 

Lemma 18. There exists such a constant G such that for each constant c, if u is a harmonic function in 
Bcr{po), with |u(p)| < cR,p G Bcr{po), then u is linear in the sphere with the radius R and the center po- 
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7.4. Relaxation process near the boundary of a slope {p,q). Now we come to one of the the bottom 
level arguments of this paper. Let us fix p, g € N. The only thing we need to demonstrate is why the sequence 
of waves near the boundary of a slope (p, q) produces a new edge of a tropical curve with the same slope, see 
Figure The set of colored vertices constituting this edge possesses the property that after waves it moves 
without changing the pattern, that is why the name self-reproducing. Once we prove this property, we can 
essentially speak about operators Gp, see Proposition All this is done via demonstrating that the toppling 
function near such a boundary is of type —qx -\- py c. 

Consider the sandpile problem on the square lattice T on the plane with a boundary ST = {(i^j)\iq > jp}. 
Let us send n waves from the point (0, N) with N very big. Let be the number of topplings at a 

vertex (i,g) during this process. We denote by the obtained state. Since the boundary is periodic we 
can model this situation by the following half-infinite cylinder. 

To be specific, we construct an auxiliary graph S, where 

y(s) = [o,p] X [0,-too) nz^ 

is its set of vertices. In this graph we connect (i,j) with (* — 1, j) for i > 0, (z, j) with (z, j — 1) for j > 0 and 
(1, j) with (p, j -I- q). Let i9S be a subset of [0,p] x [0, g], we suppose that (0, z), (j, 0) G for all z,g. 

Let us try to solve the problem of finding iJ : S —>• Z>o with AiL > 0 and H\qy. = 0,H{i,j) = n, for 
z > n -|-p -|- g. Clearly, the resultat h" of sending n waves produces the minimal solution of this problem. By 
abuse of notation we suppose that /z* are defined on E, which cases no problem since /z*(z,j) = h^{i-\-p,j + q)- 

Definition 21 . In what follows by Ca,b,... we define a function which depends only on the list a,b,... of 
variables. We use this notation for constants numerously appearing below. 

Theorem 10. There is a constant k = Cp^q^g^: G ^ such that the function n) = min(pj — qi k, n) is 

equal to almost everywhere. Namely, there is a constant Cp^q such that for the defect set 

Def = i7(i,j,rz)}, 

we have |De/| < Cp^q for each rz G N. Moreover, Def consists of two parts Def = B E where B C 
[0,p] X ]f),Cp^q\ and E C [0,p] x [a, a -b Cp^q] for some a. In fact, E represents the tracery of the discrete 
tropical edge in the direction (p, g) and B is the boundary artifacts - i.e. the colored vertices which inevitably 
appear near the boundary, but in a negligible amount. See Figure for example of discrete tropical edges 
and an example of the boundary artifacts. 

The proof of this theorem relies on two observations. The first is that the amount of colored vertices 
does not exceed a constant while rz grows. The second observation states that if a non-negative harmonic 
integer-valued function is less than a linear function on a big territory then this function is also linear. Recall 
that we call a vertex v colored if 4>{v) < 3 in the currently considered state (j). Now we derive some benefits 
of this theorem, and we postpone its proof till the end of this section. 


7.5. Dynamic and static of discrete tropical edges. Theorem [TO] says that after n waves the number 
of topplings near the boundary is equal to the function min(—gz -\-pj -\- k,n) almost everywhere. 
This provides the other point of view on the discrete tropical edges. 

Indeed, we start from the function h{i,j) = min(—gz -bpj -b k,n). Easy computations show that near the 
line —gz -bpj -b fc = zz we have Ah > I. This suggests the following procedure. We take a point (uq) where 
{Ah){vo) > 1 and construct hi which satisfies the properties hi{vo) = h{vo) — 1 and hi{v) = h[v) for all 
others points v. Then, if (A/zi)(z)i) > 1 at a point vi we produce /z 2 which is hi decreased by one at vi. 
Clearly we have the inequalities hk > h^'^'> while we find points Vk with {Ahk){vk) > 1. Theorem 10 implies 


that in the half-infinite cylinder S we have h^ = for some N = Cp^q^gs, since we need to decrease h only 
in a finite neighbourhood of the boundary and in a finite neighborhood of the line —gz pj k = n. We call 
this process smoothing of min(—gz pj -b k,n). 

Recall that a tracery is defined by the set of points where A of the toppling function is not zero. 


Proposition 11. The tracery of smoothing of min(—gz-bpj-bA:, n) and the smoothing of min(—gz-bpj-bfc, rz-bl) 
differ by parallel translation. 

Proof. These functions are differs by a change of coordinates, which is a translation (z, j) —> (* + Ci, g-bC2). □ 
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This shows that the discrete tropical edges are self-reproducing, i.e. sending a wave does not change their 
tracery. 

Thus, the discrete tropical edges emerge from smoothing of piece-wise functions which are linear on the 
one side and constant on the other. 

In general, if the toppling function is ai + bj + c on one side and pi + qj + r on the other, then we can start 
again from h = min(ai -I- hj -\- c,pi + qj + r) and smooth it. The following proposition is, then, clear. 

Proposition 12. Provided {a — p,b — q) is a primitive vector, this smoothing produce the same tracery as 
the smoothing of min((a — p)i + {b — q)j + c,r). 

Note, that total integral of is preserved, that is why the area of tropical edge in E is no more than 
(p + q)^ (see [H]). Indeed, for h = min{—qi -I- pj + k, n) the total integral of over E can be estimated 

via the sum of the values of Ah near the line —qi + pj + k = n. Recall, that we consider everything in E, so 
we will have about p + q terms of order p + q the maximal. 

7.6. Proof of Theorem llOl 

Remark 10. Denote by the amount of sand that is lost after n waves, i.e. the sand which have fallen 
into 9E. Therefore 

v^dT,wGn{v) V 

Lemma 19. The number is bounded by a constant Cp^q^QT. which does not depend on n. 

Proof. Consider the function h{i,j) = min(—qz + pj + pq,n). Clearly, it satisfies Ah > 0 and h\g-s = 
0, h{i,j) = n, for all i > n + p + q, i.e. the sand-solution h*^") is less or equal than h everywhere. Denote the 
set [l,p] X [0, 2q+\]r\ 1? by A. Since n(y) C A for ?; G dG, we have 

< 4 ^ {-qi+pj +pq) = Cp^q 

for arbitrary n. □ 

Now we consider how the picture is changing when we increase n, i.e. send more and more waves. Since 
the amount of lost sand is bounded by L, the number of colored vertices is also bounded by L. 

Proposition 13. For n big enough, there will be no topplings near the boundary and the picture will stabilize 
there. 

Proof. Suppose that the number of topplings at a vertex {i,j) tends to infinity, when n —>■ oo. Then the 
number of topplings at all neighbors also tends to infinity. Then the number of topplings near the boundary 
is also tends to infinity, that contradicts to Remark □ 

So, we can divide the colored vertices into two parts - those who stabilize and those who change their 
positions. Take the territory Tg which lies above the stabilizing colored vertices. We rely on the following 
lemma. 

Lemma 20. Let A G lA and / : A —>■ Ng be a discrete harmonic function (i.e. A/ = 0). Let f{x,y) < 
px + qy + r with p,q,r G Z. Then, li A = [0, N]^ and N is big enough (i.e. N > Cp^q^r) then / is linear on 
[Gp^q^r^ N ~ Gp^q^fj^ . 

Proof. This can be proven via direct calculation of the Poisson kernel, as in Section 2.2.3 of m- The more 
elementary approach is the following We sketch the idea and leave the details for the reader. For a harmonic 
function / we define var{f,n,A) = max(i /(z-|-n, j) — f{i,j)- Trivially, var{f,n, A) is monotone by A, 
and we prove its convexity, i.e. 2 • var{f, n, A) > var{f, rz — I, A) -|- var{f, n -|- I, A) by considering terms of 
A/ at (z, j) and (z -I- n, j). Since / is integer-valued we have var(f, I, A) > I and f{x, y) <px + qy + r implies 
that var(f,n, A) < np + 1. Therefore var(f,n,A)/n attains its maximum for some finite n = Gp^q^rj that 
finishes the proof. □ 

^We thank Misha Khristoforov for pointing us out these classical ideas. 
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Proposition 14. In the territory above the stabilized part the function is —qi + pj + c with some 
constant c. 

Proof. Indeed, in this territory the function is discrete harmonic function, and < —qi + pj + pq. 
Since this territory is growing while n grows, we can use the estimate for discrete harmonic functions: if 
is less than a linear function, then it is linear itself. Since our function is periodic and integer, we easily see 
that = fc( —qi + pj) + c with some integer k. Clearly k = 1. □ 

Proposition 15. The part with moving colored vertices is bounded, i.e. it is contained in the set [0,p] x 
[co(n), co(n) + ci], where ci does not depend on n and depends only on p and q. 

Proof. Indeed, imagine that it is possible to have a big territory above the territory where = —qi +pj + c, 
and is not everywhere constant there. We repeat the same arguments as above and see that in this case 
is a constant or of the type —qi +pj + c there. It is not possible due to the maximal principle. Namely, 
let /i" is a constant in the middle territory. Take the point where = n with a neighbor with value of 
less than n. Then this neighbor must have a neighbor with lesser value, etc. That indicates that we can not 
enter the middle territory because there the vertex have no neighbours with lesser values. We arrived to the 
contradiction. □ 

So we proved that the function come as the smoothing of the function h{i,j) = min(—gi +pj +pq, n). 
Thus we proved that an edge will detach. 

We drop a grain of sand to a point p which belongs to a discrete tropical edge. If p was coloured, then we 
have nothing to do. If not, then p has been dropped to a point with three grains and we must proceed by 
relaxation. 

Lemma 21. This relaxation has an effect only in the constant neighborhood of p which does not depend on 
the scale N. 

Proof. Indeed, remember the process of smoothing of a piecewise linear function. We start from h{x,y) = 
min(aa: by c^px + qy + r) and then start to decrease h at the points where A/i > 1. If during this process 
we decreased p, then it is easy to recover the result after the above relaxation. Indeed, let us decrease the 
function at p one less time than we do normally. Because of that we also will not decrease it at some other 
points, and that gives the result of the relaxation. It is clear, that this affects the picture only locally. If we 
didn’t change p during the smoothing, then the toppling function at p is still ax+ by + c oi px + qy -\-r what 
is not possible by the following reason. Take the territory where p belong to. Since in all neighboring to this 
territory points the function h was decreased, we also can decrease h on this territory keeping the condition 
A > 0 preserved. □ 


8. Smoothing a corner 
Consider a function /g: —)■ Z given by 

(6) fo{x,y) = min{0,pix + qiy,P 2 X + q 2 y) P Pox + q^y, 

where Pi^qi € Z and piq 2 — P 2 qi = 1- It is clear that /g is superharmonic, i.e. A/g < 0 everywhere. Denote 
by C a set of points where A/g < 0. For any r and A C denote by B^A the neighbourhood in Z^ of A of 
radius r, i.e. 

BrA = {i; G Z'^\3w G zl — raj < r}. 

The following classical theorem will be applied several times, (this should be written in the section about 
edges) 

Theorem 11. [see 0] Consider a function / which is harmonic and non-negative on a disc Bj.{0,0). Then 

1/(0, 0) - /(I, 0)1 < - niax/(?;). 

r |«|<i? 

Lemma 22. Let /: Z^ —>■ Z be a superharmonic function and there exist a number n > 0 such that the 
inequality /g — n. < / < /g holds at every point. Then there exist r depending on n such that h = /g — / is 
constant on each of the three connected components of the complement to B^C. 
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Proof. Consider a point Vg € lP‘\C such that 


Ah{vo) = -Af{vo) = 4/('Uo) - X! ^ 

vi^n{vo) 

Then there exist vi G n{vo) such that /(vo) > /(vi). If vi is not on C then again 

Ah{vi) = -Af{vi) = 4:f{vi) - /(no) - X! > 0. 

D2Gn(t;i)\{t;o} 


And therefore there exist S n{vi) such that /(vi) > f{v 2 ). This procedure defines a sequence Vk S Z^\C' 
of points such that Vk+i is a neighbor of Vk and f{vk) > f{vk+i). By adding a linear function to fg we can 
assume that fg is zero on a connected component of lP'\C that contains no- This implies that the sequence 
has at most n terms. In particular, the lattice distance from vg to C is less than n + 1, and the function h is 
harmonic outside the neighbourhood Br'(C) of radius r' > V2{n + 1). 

Now we can apply the estimate (see Theorem 111 on the derivative of h in some harmonic region. Consider 
a point V in the complement to the neighbourhood Bj-C of radius r = (v^ + l)(n + 1). Then h is harmonic 
on a neighbourhood Bn+iv. Moreover, since h is non-negative. 


\dxh{v)l \dyh(v)\ < < 1, 

where dxh and dyh are the discrete partial derivatives of h. Recall that h is integer valued and, therefore, 
constant on B^+iV. □ 


Lemma 23. Let / and g be superharmonic functions. Then h = min(f,g) is superharmonic. 
Proof. Let z; be a point. Without loss of generality, suppose that f{v) < g{v). Then 

Ah{v) = —4f{v) + E <A/(,)<0. 

wGn{v) 


□ 


Definition 22. For any n S N consider a set Pn of all superharmonic functions / such that fg — n < f < fg 
and / vanishes outside BrC for some r > 0. Define /„: —)■ Z to be 

fn(.v) = rnin f{v). 

Jfc*' n 

(the definition is much more general, can be moved to edges) 

Lemma 24. The function /„ belongs to Pn- 

Proof. By Lemmathe function /„ is well defined /„ and vanishes outside some neighborhood of C. For any 
point V there exist / S Pn such that /„(z>) = f{v). Therefore, /„ can locally be represented as a minimum of 
a finite number of superharmonic functions and so by Lemma is also superharmonic. □ 

Lemma 25. Let / and g be two integer-valued superharmonic functions such that f > g. Suppose that the 
difference h = f — g is bounded and has a maximum m. Consider a function h taking 1 on the set h~^{m) 
and 0 elsewhere. Then f = f — h is superharmonic. 


Proof. Let z; be a point where h is not maximal. Then 

4/(11) = 4/(z;) > ^ f{w) > ^ f[w). 

w^n{v) w^n{v) 

On the other hand, if h(v) = m then 

4f{v) = 4g{v)+4m = g{v) > ^ (giw) + m) > ^ f{w). 

wGn{v) wGn{v) 


□ 


Lemma 26. Consider a function hn = fn — fn+i- Then /i„(z;) < 1 for any z; G Z^. 
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Proof. It is clear that /„ > fn+i at every point of Z^.If the inequality doesn’t hold then the maximum of the 
function hn = fn — fn+i is at least 2. By Lemma 


25 


the function fn = fn — hn is superharmonic. It is clear 
that /o > /n ^ /n- If 'y is a point where h is not maximal then fn{v) = fn(v) > /o(r’) — n and if h takes its 
maximum at v then fn{v) > fn+i{v) + I > /o(ii) — (n + I) + I = fo{v) — n. Therefore, /„ S Pn and we have 
a contradiction since at some point /„ is strictly less then fn- □ 


Remark 11. It follows from the proof of Lemma 26 that the function fn+i can be characterized as the 
point-wise minimum of all superharmonic functions f < fn such that /n ~ < 1 and fn — f vanishes outside 

some neighborhood of C. 

We say that a function g is monotone in the direction e G if for any v G 

g{v + e) <g{v). 


Lemma 27. If /g is monotone in some direction then /„ is also monotone in the same direction. 

Proof. Suppose that /„_i is monotone in the direction e gZ? and /„ = fn-i — ^n-i is not monotone in this 
direction. Denote by A a set of all v G Z"^ such that fn-i{v — e) = fn-i{v), hn-i{v — e) = I and hn-i{v) = 0. 
The set A is non-empty and belongs to some neighbourhood of C because /„ is not monotone, hn is supported 
on a neighbourhood of C and hn(A — e) = {1}. Consider a function 

/; = /n - E 

It is clear that A/^(w) < Afn{w) < 0 for any w G Z'^\A. Therefore, A/^(w) > 0 for some v G A since 
otherwise we would have a contradiction with Remark II11 But 

4/„(ii-e) = 4/;(i;) < E E - e) < E Mw-e), 

w^n{v) 'w^n{v) w^n{v) 

which is imposible since /„ is superharmonic. □ 


The main result of the section is the following 
Theorem 12. The sequence of functions /„ stabilizes. 


For the rest of this section we assume that /„ doesn’t stabilize and we will show that this leads to a 
contradiction. 

Consider the functions Si{x,y) = min{0,pix + qiy),S 2 {x, y) = min{0,P 2 X + q2y) and S 3 {x,y) = inin(pix + 
Qiy,P 2 X + q 2 y), where pi and qi are the same integers as in^ Let Si be the smoothing (see ...) of Si and /g 
be the function given by 

fo{x,y) = min{s 1,82,83) +PoX + qoy- 


Clearly, /g is superharmonic (by Lemma 23) and /g 
k> N. 


A /o A /o for some A > 0. Therefore, fk < fo for 


Lemma 28. For k > N the support of hk = fk — fk+i is finite. 

Proof. If the contrary is true then there are infinitely many points where hk is non-zero. In particular, there 
is a ray of C going in the direction parallel to a primitive vector e and having infinitely many points of the 
support of hk. We can cut the neighbourhood of this ray into parallelogram pieces with sides parallel to e and 
e' where e' is such that e • = 1. Up to a translation by e there are only finitely many possible restrictions 

of the function fk to such a parallelogram. Thus, there exist two parallelograms with the same restrictions of 
fk. This allows us to define a new superharmonic function g which has a period me for some integer m and 
coincides with fk along a piece of the neighbourhood of the ray containing the two parallelograms. 

Therefore, we have dehned a new smoothing g of the edge in the direction e which is strictly smaller than 
the canonical one. This gives a contradiction. □ 

Lemma 29. There exist a point vq such that /„(ug) = /g — n for all n. 

Proof. By the assumption all the functions /„ are different. In particular, /g — /„ has a maximum n. Thus 
the supports of functions ho,..., /i„_i have a non-empty intersection for any n. Therefore, there exist vg in 
the intersection of all supports of all functions hn. Clearly, /„(fg) = /g(r'g) — n for any n. □ 
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9. Finiteness of the tropical dynamics 

In this section we will prove that we need only a finite number of steps in the recursive definition for 




(see section 4.2). 


Lemma 30. Consider an integer sequence {ki}^^ such that ki G [l,u] and for any I G [l,n] the sequence 
ki takes value I infinitely many times. Then the sequence of functions Fn = Gp^ ... Gp^ F converges to 


G. 




^F to point-wisely. 


Proof. Each term of the sequence is represented by a tropical polynomial and from Lemma the Newton 
polygon for each of this polynomials is a subset of the Newton polygon for Gp^^,,,^p^F, which is proven to 
exist by LemmaSo, in fact we need to prove the convergence for a finite collection of bounded coefficients. 
Each Gp increases a certain coefficient (see section 4.4). This implies the convergence for each sequence of 


coefficients and a point-wise convergence for a sequence of functions F^. Denote the limit by Fac,. 

Now we fix Z = 1,... ,n and take a monotone sequence rrii such that krm = I- Then the idempotency of 

for all i. Since the operator Gp, is continuous (see Remark 14) we conclude 


Gp, implies that Gp^F, 
that Gp^F^ — Too. 


= F„ 


The last equality is equivalent to the property for F^o of being not smooth at the point pi. Again using 


Lemma 30 and an implicit description for Gp 


(see section 4.2) we conclude that 


Goo — Gr, 


□ 


To proceed, we need the following intermediate result. 

Lemma 31. If pi,... have rational coordinates and F is restriction to A of a tropical polynomial with 
rational coefficients then to compute Gp.^^,,,^p^F it is required to use each Gp,, only a finite number of times. 

More precisely, any particular way of computing Gp,^_,,._p^F via Gp,, is finite. 

Proof. Eirst of all we can reduce a problem to the case when all the points and coefficient are integer. This 
can be done by a simple rescaling of A. Namely, denote by Q the product of all denominator for coordinates 
of Pi and all denominators of coefficients of F. Then we consider the points Q ■ pk on Q ■ A and a new function 
F ■. X ^ QF{Q~^x). 

Now, all objects are integer and they remain be integer while applying Gp. So, we consider GpF, it is clear 
that Gp^^,,,^p^F > GpF. It is also clear that if GpF and F are not the same, then Gp have increased some 
coefficient for F by an integer. Thus each Gp increases the sum of all coefficient at least by 1. So the process 
cannot be infinite. □ 


Theorem 13. Eor any collection of points pi,... ,p„ C A and a tropical polynomial F, the function Gp^^,,,^p^F 
can be obtained as a finite product of Gp,. 


Proof. We proceed by induction on number n of points. For n = 1 there is nothing to prove. Let us suppose 


that there is an infinite chain G„- Gr,. G„ 


which converges to G 


Pl,...,Pr 


F. Let Fi be Gp,^ Gp^^_^ ■ ■ ■ Gp,, F 
Now we take the point Pn and 
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such that Fi is very close to Gp,,,.._p„F. We can choose such m by Lemma 
stop applying Gp^. By induction we can apply Gp. to Fi with i < n only finite number times, so such a way 
we obtain F 2 . These applications can not change Fi a lot, since it is close to the limit. Now we remember 
about Pn. Let e be the minimal weighted distance from p„ to the tropical edges. Then, Gp„ applied to F 2 
increases some coefficient of F 2 by e. We claim that all the Gp, applied after that also increase coefficients 
by e, and this implies that this process can not be infinite. 

Since F 2 is close to the limit, then during the application of Gp. there will be no changes in the tropical 
curve except local ones: an edge can move slightly back and forth, a vertex can produce several vertices or 
even a new face, but in very small neighbourhood. Using the definition of Gp in Section 4.4 we observe only 
changing of the coefficients of F 2 . Let us take any pi. We can suppose that it is (0, 0) and look how we 
change coefficients. Whether (0, 0) is a vertex or an edge, all the tropical monomials around it are of the type 
ax + by, without a constant; then by the same reasoning as in Lemma we see that all incrementation are 
multiples of e. 
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Therefore, starting from Fi we applied only finite number of operations and arrived to the limit. That 
means, taking N big enough that Gp.^ i ''' +i contains all these operations, we also arrive to the 
limit. □ 

Remark 12. Theorem [T^ implies Lemma 


10. One point case 

Now we are going to discuss in more details what happens if we add one point to the maximal state. The 
relaxation for the state of this type serves as a local model for more general relaxations. In particular it 
provides an easy description for the operator Gp. Another goal of this section is to show how to apply in the 
most basic situation the general techniques we developed. 

For the case when the underlying graph of the model was a square piece of a rectangular lattice the result 
has been already described in the first section (see Figure [^. Here we continue to work with a subgraph of 
the same lattice but the boundary now can be a lattice polygon of an arbitrary shape. 

Lets fix 5, a large enough polygon with rational slopes. Consider another polygon 5' which is constructed 
as a union of 5 with its two horizontal and two vertical shifts 

5 u 5 + (±1,0) U (5 + (0, ±1). 

This new polygon can be seen as a set of those points in that can be moved to <5 by the action of standard 
generators of C (or already were in <5). Now we define a graph F = r(5). 

Definition 23. Let r(^) be a full subgraph of 1? whose set of vertices is 5' . The boundary of this graph 

is defined as (^'\(5) n I?. 

Consider the maximal state (j)Q and a vertex v € V(r)\dr. Let us add one extra grain of sand to (f>o at 
V obtaining an unstable state </>. Suppose that v is far enough from the boundary of S. This quite vague 
assumption together with the assumption that <5 is large enough could be made more rigorous in the following 
way. In the spirit of the second section we fix a lattice polygon A and a point p in its interior. Then we can 
take a positive integer N and define <5 to be A • A and v to be the coordinate-wise integer part for N -p. Then 
as we claimed before in the Theorem [l] we expect to see a very regular behavior for the relaxation of (j) when 
N is large enough. 


10.1. Corners of A are smooth. Now we suppose for simplicity that all the corners of A are smooth, i.e. 
the primitive vectors for the adjacent edges span a parallelogram of area 1 (see definition 241. The general 
situation will be described in the end of the section. 


Definition 24. A corner of a polygon is called smooth if 

det\wi,W2\ = ± 1 , 

where wi and W 2 denote the the primitive vectors in the directions of edges adjacent at the angle and \wi,W 2 \ 
denotes a 2 x 2 matrix made of the vectors wi and W 2 - A polygon is called smooth if all of its corners are 
smooth. This terminology is inspired by the theory of toric surfaces. 


To understand the relaxation for (ft we decompose topplings on waves. As it was explained in the Proposition 
13 after applying a number of waves the picture will stabilize near the boundary of the polygon 5. At his point 


the set of colored vertices will consist of boundary artifacts, edges parallel to the boundary and edges going 
from the corners of the boundary. Each new wave emitting from v doesn’t touch boundary artifacts, shorten 
edges parallel to the boundary and move them towards v (see Lemma 32) and prolongs edges going from the 
corners. During this process some of the edges can collapse. This creates a vertex outside the territory of v. 
The collision is shown in the Figure]^ for v = [A • p]. 


Lemma 32. A periodic pattern in the primitive direction (a, fi) after applying M waves moves at the distance 
M{a‘^ + i'^)~^G _|_ o(A) towards the source of a wave. 


Proof. Without loss of generality, suppose that a > 0 and /3 > 0 and the pattern is arranged along a line 
passing through (a, /3) and the origin. Lets send M waves from the point (A, 0). Let h be as usual the height 
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function for this process. Then in the territory of the source of wave it is equal to M. Thus the new location 
for the pattern is along the locus of non-smoothness of the function 


{x, y) I—)■ min{M, fDx — ay). 

This locus is again a line with the same slope. The new line passes through a point (M//3,0). And so the 
distance between the lines is 


(n,(M//3,0)) 


M 

+ p2 ’ 


where n is an orthonormal direction for the pattern. 


□ 


We continue sending waves until one of the edges parallel to the boundary will not pass through v. This 
simply means that at this moment v went outside its former territory and waves can no longer be applied. 
To complete the description of the relaxation for ^ we simply add one grain of sand at v to our current state 
(we can do it safely because v doesn’t belong to any territory). From this description it is clear that the 
dependence on JV is quite weak and we can actually find a limit for the picture when JV tends to inhnity. To 
do this explicitly we introduce a weighted distance function. 



Figure 10. A result of sanding the same number of waves from a point v near a non-smooth 
corner and its perturbation. The non-smooth corner is spanned by the vectors (—1,2) and 
(2, —1). On the second picture we cut the corner with an edge in the direction (1, —1). 


Definition 25. A weighted distance function Za is a function on points p in a lattice polygon A. It is defined 
as a minimum of dist{p, e)\lo? over all edges e in the boundary of A, where (a, /3) is a primitive vector 
in the direction of e and dist(p, e) is a distance from p to a line containing e. 

The function /a measures how fast a point will meet the set of colored vertices. It also appears that it 
counts the number of topplings at v. 

Proposition 16. Let F)v be the toppling counting function during the relaxation of the maximal state on 
r(N ■ A) after we add an extra gain ai v = [N ■ p]. Then 

^Fn{[N -p]) —> mm{lA{p),lAip)) 

10.2. General A. When we apply waves near a non-smooth corner it produces more edges then in a smooth 
case. An enlightening way to treat a non-smooth corner is to consider it as a collection of degenerate smooth 
ones (see Figure [To|). 

Except for this pathology of extra edges coming from bad corners, nothing have really changed. Concep¬ 
tually, we can thing about a general polygon as a degeneration of a good one, this principle is realized in 
Proposition After that, we start to emanate waves from the point v as we did before. And we do that 
until some of the edges parallel to an edge of the perturbation of the dA passes through v. This procedure 
dictates us the general definition for the weighted distance function. 

The function Zq serves as an example of a tropical power series. In this terms the interior of the convex 
hull for n should be thought as a domain of convergence for Iq. We will discuss more of that in Section [T^ 
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With the current definition of Ia Proposition!^ works for a general A without any changes. What is more 
questionable is a consistence of Definitions and and potential infinitude for I a ■ This two problems are 
resolved by the following lemma. 

Lemma 33. Let P(A) be a set of all primitive vectors in the union of all triangles with vertices n(ei), n(e 2 ) 
and 0, where ei and 62 are adjacent edges of A. Then 

ln{x) = inf {w-x-c^) 

Here n(e) denotes a primitive normal vector to e towards the interior of A. It is clear that T’(A) is finite. 
And for A having only smooth corners T’(A) coincides with a collection of all primitive normals to the edges 
of A, i.e. P(A) = n{E(A)). 

Proof. It is clear that non-primitive vectors do not contribute to I a because kw ■ x — Ckw = k{w ■ x — Cw) and 
w ■ X — Cyj is nonnegative for all x € A. 

Let ui be a primitive vector not orthogonal to any of the boundary edges. This implies that tc • x — c^, 
vanishes exactly at one vertex of A. We can move the origin to this vertex (which makes c^, to be equal 0) 
and consider the two edges ei and 62 going from the vertex. So we have w = ain(ei) -I- 0211 ( 62 ). And ai are 
positive because ai\n(ei)\'^ = w ■ Ci > 0. 

Suppose that oi -I- 02 > 1, i.e. w is outside of a triangle with vertices 0, n(ei) and 71 .( 62 ). We want to prove 
that w will not contribute to I a- The whole polygon A is a subset of the cone spanned by ei and 62 . So it 
is enough to show that w- dominates either ?i(6i)- or 71 ( 62 )- in the cone. Let p be a point in this cone, i.e. 
71 ( 61 )-p and 71 ( 62 )-p are nonnegative. Then w-p = ai7i(6i) •p-|-a27i(e2)-p > (oi-|-a2)77iax(7i(ei)-p, 11 ( 62 )-p) > 
max(7i(ei) • p, 71 ( 62 ) • p). D 


11. Sand approximation for tropical polynomials 


In the next two sections we are going to show how statements of Section]^ can be deduced from the results 
of Section For any tropical curve we will construct its sand model called a discrete approximation. This 
model will be represented by a particularly nice set of colored vertices on a large lattice. In this approach the 
operators will be interpreted in terms of certain relaxations process. In the previous section we saw 

the most basic realization of this principle for the case of an empty curve and Gp. 

Let A be a lattice polygon and F’ be a non-constant tropical polynomial vanishing at the boundary of A. 
Let C be the intersection of a curve given by F with A. Consider a graph T = T{N ■ A) (see Definition 10) 
for N large enough. 

We want to construct a new state </> such that its set of colored vertices, i.e. where </> < 4, contains N -C in 
its tubular neighborhood of radius m, where ttiN is fixed for given F and doesn’t depend on N. If we have such 
a state, then any connected component in the complement to a curve C would naturally define a territory 
for (j). As another requirement for (j) we want this mapping to be injective, i.e. the territories separated by 
the edges oi N ■ C should be distinct. The territories for (p which are in the correspondence with connected 
components of A\C will be called large territories. 

Last requirement is necessary to leave the state untouched outside a small neighborhood of a large territory 
while dropping a grain of sand in the territory and then relaxing. During this process we could have certain 
artifacts near the boundary of the territory. Clearly, there are many states p satisfying the conditions listed 
above. And if we chose a generic p the artifacts will inevitably appear. 

To avoid this problem we could specifically choose the structure of p along an edge e in C. Namely, we can 


fill a neigborhoud of A • e with a periodic pattern corresponding to the slope of e (see Section 7.51. This will 
ensure that the colored set along N ■ e moves without artifacts during the relaxations. And this is almost it. 

It only remains to take care of the vertices. To avoided artifacts near trivalent vertex we can model (p by 
a corner patch arising near the corners of the boundary. In the case when our curve has vertices of higher 
valence or cycles we have two issues. 

The first issue can be easily resolved. We can consider a tropical curve N ■ C (or rather a function 
p I—>■ F{j^p) for p £ N ■ A) and then make it generic using a perturbation of order -s/N. Thus all the vertices 
will be trivalent and macroscopically picture remains the same. 

The second issue is more delicate. We need to show that we can actually glue all the patterns for edges 
with the corner patches while going around the cycles. Indeed, we will see that this is possible. The idea of 
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the proof is very simple. Instead of making a surgery on patterns we can naturally grow a whole curve by 
emitting many waves from different places. We formulate it in the following proposition. 


Proposition 17. For any N there exist a smooth lattice polygon S^, such that all are dual to the same 
fan S{Fan{A)) (see definition 26), length for each edge tends to infinity and -^Sn converges to 6, and a 
sequence of stable states 


on r((5jv), where each 4>i+i is a result of sending a wave for (j)i and a sequence of toppling number functions 
/iTV: —>• Z of the processes (j)o —t converges to F, i.e. 


h^{[N.]) F(.). 

iV Af—>00 


The state (j)o is as usual a maximal stable state and as before counts the number of topplings at 
each point while transforming (j)Q to 0™'^, outside V{T{Sn')) the toppling number function is naturally 
prolonged by 0. The proposition can be seen as certain inversion for the Theorem]^ 

In the statement of the proposition we used a notion of a fan and a concept of duality for fans and polygons 
that were not introduced yet in this paper, here is a definition for them. 


Definition 26. A (complete, lattice) fan $ is a finite collection of rays with rational slopes starting at the 
origin such that the angle between two consecutive rays is less than tt. To encode a fan we take a set of all 
primitive vectors in the directions of its rays and say that a fan is generated by this vectors. A polygon A 
is dual to a fan $ if $ is generated by all the vectors n(e) for e G E{A), where n(e) is a primitive normal 
to vector to an edge e (see the page for the definition of n(e)). Let S{^) denote the (minimal) smooth 
refinement for iS'($) that is a fan generated by all the primitive vectors from the triangles spanned by 0, wi 
and W 2 , where wi and W 2 are the primitive vectors in the directions of two consecutive rays in $ (see figure 

TTl) 



Figure II. Two polygons and their dual fans. The fans are generated by collections of 
primitive vectors, that are shown on the picture. The right fan is a minimal smooth refinement 
of the left one, and the right polygon is a result of corner cuts for the left polygon. 


Proof. Fix N ^ 0. First we need to define a polygon Sm- Consider a lattice polygon [-s/iVjA. We can make 
each of its vertices smooth (see Definition 24) by doing a sufficient number of small lattice corner cuts (see 
Figure 10). This gives us a polygon A)y dual to the fan S{Fan{A)). We define a polygon Sn as [V^JA'^ 
and graph Fjv as r(5Ar, 1). Clearly jf6N converges to A in the Hausdorff sense. 

Consider the state (fo- We start deforming it by sending [-x/iV] waves from the middle of A • A. This will 
create one territory of the area proportional to in the center of the polygon. Also for each edge e of A 
we will have a territory with area proportional to A^/^ adjacent to (the neighbourhood) of A • e. And finally, 
for each non-smooth vertex u of A we will see a number of territories with area of order 0{N) near A • v. All 
other territories will be of the area 0(1). 
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Consider an edge e of dA. We know that F is linear in the neighbourhood of e and vanishes at e. This 
means that near e the function F can be written as 


F(x, y) = m{e){qx — py c) 
for m[e) € Z, coprime p and q and c G K. 

Let Cl ... Cm be the sequence of all edges in A. Take a territory near N ■ ei with m(ei) > 1. This territory 
has a shape of trapezoid. More formally we can say that the territory sits in the small neigbourhood (for 
example of radius of the trapezoid. The larger base of this trapezoid is A • ei. Consider a point in the 

territory near a smaller base of the trapezoid. We emit ['/N/m{ei)\ waves from this point. This will split 
the territory on two territories with areas of order Both of this territories are again trapezoids. 

If m(ei) > 2 we repeat the procedure for the trapezoid adjacent to the boundary. This will again split the 
trapezoid into two pieces. Thus we create m(ei) — 1 consecutive trapezoids with bases parallel to ei. 

After we pass to iV • 62 . We do exactly the same thing but the number of waves emitting from one point 
is equal to [ m(e|i^(e 2 ) ]' general when creating a trapezoid territory for Ck we use m{ek) times less waves 
then for a single territory parallel to Ck-i- 

At this point all the colored vertices are arranged along the boundary and a small neighbourhood of a 
certain tropical curve Cn- Note that the territories with the area 0{N°‘) for a > 0 correspond to the connected 
component of the complement to Cn- 

The main benefit of our construction oi 6^ is that is smooth. This implies that the connected compo¬ 
nents of the complement to Cat, are in one to one correspondence with the integer points inside its Newton 
polygon. This means that by sending waves we can freely deform all the ’’coefficients” for h^. Since the 
Newton polygon for Cn contains the Newton polygon for our initial function F we can deform Hn to those 
that converge to F. □ 


Definition 27. A sequence of polygons 6n, functions Hn and states <j)N 
Proposition 17 will be called a discrete approximation for F. 


satisfying the conditions of 


It is clear that determines (pN and F. So abusing the notation we will say that is a discrete 
approximation for F. 


Remark 13. It is clear that a choice of exponents for N in the proof of Proposition [T^ is not crucial. In 
particular, the size of the cuts for iV • A and the order of perturbation for a discrete curve could be taken 
smaller than '/N . 


12. Sand interpretation of Gp 

Finally, we are going to debunk the mystery of tropical relaxation. In particular we will show that the 
curves and functions arising from the sandpile model enjoy certain extreme properties. 

Lets fix a discrete approximation /iat for F. Take a point p G A. For each N we modify (j)jq in the following 
way. If [N ■ p\ is in the territory for (pjq we add a grain of sand at the vertex [N ■ p], then relax and remove a 
grain of sand at [A • p], otherwise we leave (pN unchanged. This gives a new state We define a function 
on vertices as a the number of topplings during the relaxation. 


Proposition 18. There exist a tropical polynomial F' vanishing at the boundary of A such that a sequence 
of functions = hpf + hjq and states (p'j^ is a discrete approximation for F'. Moreover, 

• F'{x) = Gp{F){x) for all x G A; 

• if p belongs to the connected component n of the complement to a curve defined by F then 


F'ix) 


F{x) + Iq{x) X G ft 
F{x) X G A\fl. 


Recall that % denotes the weighted distance function (see Definition 34). Here D of course is a polygon, 


so we can use lemma 33 to see that Iq is a tropical polynomial. 


Proof. The state is a result of applying waves to 4>o because we can decompose the relaxation in the 
construction of into a sequence of waves applied to 4 >n- Thus the only thing we need to prove is that 

converges. 
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Consider a vertex = [N ■ p]. If p is on a curve C defined by F then from the convergence of the 
distance between vn and the boundary of the boundary of the territory is o{N). This implies that — Hn 
is o{N). And thus h'j^ converges to F' = F. 

If p is not on C then vn is inside of a large territory T. The state is a result of applying a maximal 
number of waves at v. Each wave will simply move each edge in the boundary of the territory towards v, 
as it was described in the previous section. Each ed ge e ' of the territory T will be moved on the distance 
Mj after M consecutive waves (see Lemma 321. Here (q;,/ 3) is a primitive vector parallel to e and 
M > 0 but M < IV. 

We can actually calculate h' up to o(n) at every vertex v' . If v' is not at the same territory T then 
h'(v') = h(v). Otherwise, each wave emitted from v will contribute I to h'(v') — h{v') until v' or v is not 
at the boundary of the territory. This means that if l(v) > l{v'), i.e. v' reaches the boundary first, then 
k'niv') = h]y{v') + l{v') + o{N). Here I = In-A' denotes a weighted distance to the boundary of • A' 
(see Definition 25), where A' is a connected component to the complement of C corresponding to the large 


territory T. And if l[v) < l{v') then h'{v') = h(v') + l(y) + Of^N). It is clear that after rescaling the function 
I vanishes at the boundary of the territory T and after rescaling converges to a piecewise linear function on 
the connected component to the complement of C. This implies the convergence for h'j^. □ 


Remark 14. Proposition 18 implies that Gp is a continuous operator on the space of tropical polynomials. 


13. Convergence for the sand relaxation 

In previous sections we developed some kind of theory that provides a nice understanding of certain 
relaxation processes. Unfortunately, we didn’t have an occasion to prove several first theorems. We’re going 
to fix this. 

Consider a tropical polynomial F vanishing at the boundary of a lattice polygon A. We fix a collection of 
points pi,... ,pn in A°. 


Consider a discrete approximation for F (see Definition 27). Consider the sequence of stable states (pN 
defined hy Hn on graphs S^. We construct a new sequence of possibly non-stables states (j)% by adding extra 
grains of sand to pN at the points [N ■ p^], i.e. 


( 7 ) 




(I)n{v) + 1 V G {[Npi],...,[Npn]} 
otherwise. 


As we did before we count the number of topplings at each point v G V{Sn) during the relaxation for 
This gives a function hff : V{5pf) -G Z>o. 

Proposition 19. The sequence /ijv converges to a function on A, i.e. 


( 8 ) 




Moreover G 




F = F + F. 


Proof. Eollowing the recursive definition for of Section]^ suppose that the curve given by F passes 

through all the points pi. It means that for N large enough each of the points [N ■ pi] is in the small 
neighborhood for same edge in the set of colored vertices. So from the Lemma 21 it follows that hisf = o{N) 
and so F exists and is identically zero. 

Suppose now that some point pk is not on the curve defined by F. This means that for a large N the vertex 
[N-pk] is in a large territory of the state pN. We modify the state pN by adding an extra grain at [N-pk] and 
then relax it and remove a grain at [N ■ p^] . Denote the resulting state by Proposition says that a 

sequence of states p'jy together with corresponding toppling counting function is a discrete approximation 
for a function F'. And from Proposition it follows that F' = Gp^F. 

In the last paragraph we mimic the induction step (we will call it the step f or t he point pk) in the recursive 
On each step we contribute h'j^ to hjq. From Theorem 


definition of Go 


13 


we know that for all A 0 


the number of steps is finite, doesn’t depend on N and the steps can be chosen coherently for different N 


then we will prove and G. 


■'pi. 


F = FFF. 
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On the other hand, each step provides a partial relaxation for the state defined by Q. Since we know 


"N _ 

that any relaxation on a graph Sn should terminate, to complete an honest relaxation of 0^ for any given 
N it is enough to do only finite number of the steps and in the end add one grain at each vertex [Npi]. This 
observation also implies that the steps can be applied in an arbitrary order. 

Now we want to analyze how each inductive step for a point pk terminates on the level of discrete approxi¬ 
mations. Recall that the step consist of applying a series of waves at the same large territory where the point 
[Npk] sits. We stop sending waves when [Npk] appears to be at the boundary of the territory. 

On the limiting level this correspond to passing of the curve given by Gp^F through a point [Npk]. 


Moreover, by calculating the speed of each edge (see Lemma 32) under the action of waves we know exactly 


which edge (or vertex) of the curve will pass through the point, i.e. we index termination of steps with such 
combinatorial events. 

Since the set of colored points is located along the curve, for a large N the cause of termination of the 
step is the same as for all N. This implies that for Ni,N 2 ^ 0 after applying a sequence of steps for the 
points pki, ■ ■ ■ ,Pkm G {pi ■.. ,Pn} to the states and ^^2 the combinatorial distributions for the rescaling 
of points Pi... ,pn on big territories of c/iati and (pN^ coincide. Thus, the steps can be applied coherently to 
different cpN ior N large enough. □ 

If we apply the previous proof to the case when F = 0 we get essentially a proof for Theorems and 

Proof of Theorem^ As we explained in the introduction, Theoremis an easy consequence of Theorem 
So we will concentrate on Theorem [21 

The only difference between the statement of Theorem |2] and the essence of the proof of Lemma is that 
the theorem was formulated without the involvement of discrete approximations. Comparing with Section [^ 
in this theorem and in the whole introduction we always start with an empty curve on a lattice polygon A 
and an identically zero polynomial function. 


In this case the discrete approximation (see Definition 271 constructed in Proposition 17 is very easy to 
describe. The sequence of states (pN consists of maximal states and the toppling number function /ia? is 
identically zero. The only nontrivial thing here is the choice of a sequence of polygons Sn for the underlying 
graphs of sandpile models. 

These polygons 5n are not very different from the rescaling of A, they are essentially derived from N ■ 6n 


by performing lattice corner cuts for for non-smooth corners. As it is claimed in the Remark 13 we can take 
smaller exponents for the cuts then in the original proof of Proposition |17| 

So we construct ^ai by taking the corner cuts of order for the polygon N ■ A. Together with a sequence 
of trivial toppling counting functions hpf = 0 this defines a discrete approximation for F = 0. Applying the 
proof for Lemma to this discrete approximation we get the convergence for the relaxation of a sequence of 
states (p'j.^ (see equation 0). where i^ai is a sequence of maximal states. 

Now we need to estimate the difference between the relaxations for (p'j^ on Sn and the analogous state on 
N ■ A. Clearly <5 at is a subgraph for TV • A, so we can think of a relaxation on Sn as an incomplete relaxation 
on TV • A. Other way round, to have a relaxation on we can run a relaxation on iV • A and freeze the 
vertices from N ■ A\dN, i.e. we forbid to do topplings at these vertices. 

After the partial relaxation there can be more than 3 units of grain at each of the frozen vertices. If we 
unfreeze one grain of sand at a vertex v it will produce a certain number of waves until the amount of sand at 
V will be reduced. Each wave changes the toppling number function at each vertex of TV • A at most by I. A 
number of waves that will be emitted from v is of order since the distance from v to the boundary 

of A • A is at most So the total distortion for the toppling number function coming from unfreezing 

the grain at v has an order 0(A^/^). 

Lemma 34. There is a constant Cn^A < oo depending only on the number of initial points pi,... ,pn S A° 
and the polygon A, that bounds an amount of sand at each frozen vertex v € (A • A)\Sn after the relaxation 
on Sn. 

Proof. The existence of the constant follows from the fact that there is only finite number of different types 
of V. We will define a type of a vertex in -APSn in such a way that vertices of the same type will have the 
same amount of sand in them after the partial relaxation. 
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The simplest and probably the largest type consists of the vertices not adjacent to N - A, such vertices 
need not be relaxed. Some vertices are located around the corners, there number doesn’t depend on N ^ 0, 
so we can declare each each such a vertex to be of its own type. 

The last type consists of those vertices that are located on the cuts, i.e. the small edges ot Sn- There 
number tends to infinity with N, but there is only a finite number of locally different vertices. In this case 
we declare two vertices Vi and V 2 to be of the same type if they both are on the same edge e of and 
vi — V 2 = mw for some m S Z, where w is a primitive vector in the direction of e. □ 


The total number of frozen points has an order And from the computation of the 

distance to the boundary we deduced that each extra grain at these vertices produces at most By 

Lemma 34 a number of extra grains at any fixed frozen vertex is at most Cn,A- Thus at every vertex of the 
graph on A • A the difference between the difference between toppling number functions for the relaxation on 
5n and A • A is at most Therefore, the two toppling number functions converge to the same limit 

in the sense of equation ®. □ 


Despite the fact that the logic is a bit twisted and dependence of claims slightly counterintuitive, the rest 
of the proofs will go much easier. For example, what we’ve done during the proof of Theorem clearly gives 
the proof for Remark]^ that says that Gpi,..,p„(0) coincides with a piecewise linear function from Theorem]^ 
Lemma follows from Lemma Remark and Proposition in the following way. We consider a 
descrete approximation for the function F(a:i,..., a;„). From Proposition 19 it follows that if we 

add a gratin at to (pN and relax it then the resulting state will approximate F(xi,... ,a;„+i). Since 

the curve given by F{xi ,..., Xn) passes throug Xn+i we know that the vertex is at the distance o{N) 

from some descrete tropical edge in We what to show that it is required to send at most o{N) waves from 
[Npn+i] to complete the relaxation. This would imply that the resulting descrte approxiamtion has the same 


limit as the original one. The case when [Np^+i] is inside a descrete tropical edge is covered by Lemma 21 


So we can assume that [Npn+i] outside a tropical edge. Thus, by Lemma 32 we can estimate the number of 
waves as a constant depending on an edge times o{N). 

Theorem 1^ that claims the minimality of Gp^^,,,^p^F, actually follows from the correctness of the implicit 
description for Gp^^,,,^p^F. Together with Theorem]^ and Lemmaj^this implies Theorem]^ Lemma|^is an 
easy consequence of Lemma |21[ And finally, in this presentation Proposition is a trivial consequence of a 
definition and correctness for and Remark]^ 

Here we prove that the limiting procedure for sandpiles defined on an convex body gives almost the same 
answer as for polygonal boundary. Namely, the critical locus tends to a tropical curve, but it has infinite 
number of vertices now. For that we develop some polytope theory and study different convergency questions. 


14. Tropical series 

Definition 28. We say that a polygon A C is a Q-polygon if all its sides have rational slopes. 

Definition 29. [12] A function F : D — > K is called a tropical series if for each x G ft° there exists a 
neighbourhood U C D of a; such that F\ij is a tropical polynomial. 

Let D be a convex domain and p = {pi, ... ,p„} be a collection of points in D°. We always assume that 
the collection of points p is non-empty. Denote by V{fl,p) the set of tropical series F on D such that F|n > 0 
and F is not smooth at each of the points Pi. 

Definition 30. For a; S define Gp(x) = inf{F(a;)|F G F(D,p)}. 

Lemma 35. The function Gp is a tropical series. 

Proof. Let g be a point in . If in a neighbourhood of q the function Gp is the minimum of a finite number 
of tropical polynomials then it is itself a tropical polynomial. It remains to show that near q the function Gp 
can not be a minimum of infinite number of series non-trivialy contributing to Gp. If it is not the case, then, 
using the fact that all our linear functions are of the type ix + jy + aij,i, z S Z, we obtain that Gp{q') = —oo 
for some q' arbitrary close to q what is impossible since Gp >0. □ 

In the notation of Section l^the function Gp is Gpj_..._p^(0). 
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Figure 12. Additional grain of sand thrown to the center. Sand picture and limiting example. 


Remark 15. Each tropical series / is a concave, piece-wise linear real-valued function on 12°, in general / 
cannot be extended to ft. 

Corollary 2. The functional F Fdxdy defined on the space V[fl,p) has a unique minimum Gp. 
Definition 31. Let be a tropical analytic curve defined as the set of points where is not smooth. 


Note that Gg is a tropical polynomial if 12 is a lattice polygon, as we proved in Lemma 


33 


Definition 32. A tropical analytical curve in 12 is the locus of non-linearity of a tropical series on 12. 


Definition 33. For any w denote by Cw the infimum oi w ■ x over x € 12. Let Aq be the set of all such 
w that Cm 7^ —oo. If re S Aq define function Iq to be 

= w ■ X — Cu). 

Note that is positive on 12° and it is actually a support function for 12. 

Definition 34. The weighted distance function la on 12 is defined by 

laix) = inf l^{x) = “ V Cmx“ ” 

wGAn ' 

wGA[i 

Since la is an infimum of non-negative functions, it s a non-negative concave continous function on 12. For 
each point x S 12° 

Definition 35. Denote by G^ the tropical analytic curve defined by la- 

Definition 36. Let Ra be a set of points on x G 912 such that /q(x) = 0 for some w G Aa- 

Lemma 36. The set dft\Ra is a disjoint union of intervals with irrational slope. 

Proof. Consider a connected component / of dfl\R. Note that if / is not a straight interval, then we can find 
a line {la = 0} such that its intersection with I is not empty. Thus, / is a straight interval that cannot have 
a rational slope because otherwise I C R. □ 

Lemma 37. For a smooth convex function / with /(O) = 0 we have /((/')“^(x)) < x for x close to 0. 

Proof. Indeed, denote y = (f')~^{x), then we have f{y) < f'{y) which easy follows from convexity for 
y <1. □ 

Lemma 38. For each convex function / such that / > 0 on (0,£] there exists a smooth convex function 
/ > 0 on (0,e'],e' > 0 such that / < / on [0, e']. 
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Proof. We can consider the function fi{x) = f{x) ■ x'^, then choose a sequence —)■ 0 and then draw piece- 
wise linear /2 through points fi{pi) and smoothen it. Namely, the fact that the interval between (pi, fi{pi)) 
and {pi-i, fi{pi-i)) is below the graph of / follows from > f{pi-i)pf_i, so we choose pi = e and pi 

from the condition = 2/(pi_i)pi_i. Clearly pi < and the only limiting point p* of the sequence of 

Pi is p* = 0 since ^ = 2f{p*)p* implies f{p*) = 0 or p* = 1/^/2 but we chose e' small enough. Then we 
smooth this piecewise linear function. 

□ 

Lemma 39. The function la can be continuously extended to dVl by zero. 

Proof. It is clear that Iq is zero on the closure of zero sets of functions w € It means that it is enough 
to consider the situation near irrational smooth edges. Consider an edge of a boundary with irratioanl 
slope a. By the previous lemma we can suppose that the boundary near the edge is given by a graph of 
a smooth function ax -I- f(x). By continued fractions for an irrational number a we get two sequences of 
numbers, P 2 nlq 2 n < a < P 2 n+i/q 2 n+i such that a-p 2 n/q 2 n < l/g|„,P2n+i/<72n+i - a < l/gi„+i- Take an 
approximation p/q = Pi/qi for a. 

Then for a tangency of slopep/g we have a + f'{xQ) = p/q, then the equation of the line is y—axQ — f{xo) = 
p/q{x — xq) and we want its value at {—t, —at). Id est, xq = {f')~^{p/q — ce), 

q{y - a{f)-^{p/q -a)- f{{f)~^{p/q - a))) = p{x - {f)~^{p/q - a)). 

7 

-t{aq -p) + {p- aq){{f')~^(jp/q - a)) - qf{{f)~'^{p/q - a)) 

Since \aq — p\ < 1/g we have | — t{aq — p)\ < t/q, 

{p- aq){{f')-^{p/q- a)) < l/q{f)-^{l/q^), 
which tends to zero as g —)■ oo, and since fiif')~^ix)) < x for x small enough, we have 

qf{{f')~^ip/q - a)) < qip/q -a) < i/q. 

□ 

Remark 16. If f{x, y) = ix + jy + aij,i,j € Z, Oij € M, /|n > 0 and {i,j) (0,0), then f > la on Q. 

Remark 17. For a point g € we have = min{?Q(x), lo(( 7 )}. 

Proposition 20. It follows from the previous remark that there exists e > 0 such that Gp(x) > ln(x), for 
any x in the e-neighbour hood of dit. 

Proposition 21. For any x G il, and p = {pi,... ,Pn} 

Gpix) < < n ■ la(x) 

15. Limits of toppling functions for an arbitrary boundary 

Lemma 40. Consider a function F on V, given by F{x) = inf G^{x), where the infimum is taken over all the 
Q-polygons A containing ft. Then F coincides with Gp. 

Proof. Consider an arbitrary point q G V.°. Since G^|a € V{fl,p) for any A D O, we have F{q) > Gp{q). 
On the other hand, there exist only a finite number of monomials in Gp which contribute to the values of 
Gp at all the points Pi G p and q. Thus, eliminating all the other monomials from Gp, we get a new tropical 
polynomial Fq. Note that Fq is not smooth at the points pi and Fo(q) = Gp(q). Moreover, the set of points 
where Fq > 0 contains a Q-polygon Aq D fl, because Fo|n > 0. This implies that Fg G V(Ao,p) from which 
it immidiatly follows that 

G^(x)=Fo(x)>G^°(x)>F(x). 

□ 

Lemma 41. For any e > 0 the set fig = {x G fI|Gp > e} is a Q-polygon and Gp|o^ is a tropical polynomial. 
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Proof. Since Gp is continous and vanishes at the boundary of ft, the set Gp = e is a curve disjoint from dV,. 
We claim that the intersection of with the tropical analytic curve dehned by the series G^ is a finite part 
of a tropical curve. Indeed, if the intersection is supported on an inhnite graph, then we can find a sequence 
of vertices convergent to a point y £ V.°. Thus, there is no neighbourhood of y where the series Gp can be 
represented by a tropical polynomial, this gives a contradiction. 

The finiteness of the number of vertices implies that there is only a hnite number of monomials participating 
in the restriction of G^ to the domain □ 


Remark 18. For a convex domain we constructed a canonical series of polygon fig —>■ ft. This could lead 
to new toric insights. 

Lemma 42. In the hypothesis of the previous lemma, for e > 0 such that Gp(pi) > e for each pi £ p we 
have Gp' = Gp - e. 

Proof. We can canonically extend Gp' + e to fl, using the fact that it is a tropical polynomial. Since Gp' + e 
is linear near each side of dflg and Gp^ + £ < Gp inside fig, we have Gp‘ + e > Gp on fl \ fig. Then, 
Gp' + e S V{fl,p), that finishes the proof. □ 


Lemma 43. Consider a function F on fl given by F{x) = supGp(x), where the supremum is taken over the 
set of all lattice polygons A C fl containing x. Then F coincides with Gp. 


Proof. It is clear that Gp > F on fl. Consider a point q G fl° and e > 0 such that Gp(q) > e and Gp{pi) > e 


for each pi e p. We are going to prove that F{q) + e > Gp(g). Lemma 41 implies that the set fig = {x S 


fI|Gp(x) > e} is a Q-polygon. Evidently, q € fig and F S V{flg,p), therefore F{q) > Gp'(g) = Gp(g) — £. 


□ 


proof of Theorem We are going to prove that a sequence of functions Fp: fl —>• K given by 


FN{x,y) = -^Toppl^^{[Nx], [Ay]) 

uniformly converges to Gp. It is clear that F 7 v(A)(x,y) < Fjsr{fl){x,y) < F 7 v(A')(x,y) for polygons A, A' 
such that A C n C A'. Then, using Theorems from the previous article we get Gp = limF 7 v(A)(x, y) < 


lim F 7 v(A')(x, y) = G^ . Then, applying Lemmata 


get 
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and 


43 


by standard squeeze convergency arguments we 


limGp = limFNifl){x,y) 


limGp' =Gp. 

A' P P 


Finally we see that, the sequence of sets has a limit G in fl in the Hausdorff sense. Indeed, -ktEn- 


tends to Gp' on fig for each e > 0. 


□ 


Corollary 3. We can do all the same for all convex domains, possibly unbounded. 

Corollary 4. Furthermore, under the above hypothesis, G = C\dfl is the zero set of an inhnite tropical 
series Gp. Moreover, G passes through the points pi,... ,pn and dA \ G is the collection of intervals with 
rational slopes. 


Remark 19. Suppose that dfl contains a line with irrational slope. In this only case there is no convergency 
of the toppling function. 


16. Tropical symplectic length of curves defined by tropical series 

Definition 37 (See [20]). The tropical symplectic area of a finite segment I with a rational slope is given by 
Area{l) = Length(l) ■ Length{v), where Length^—) denotes a Euclidean length and u is a primitive integer 
vector parallel to L If G' is a finite part of a tropical curve, then its symplectic area is the weighted sum of 
areas for its edges, i.e. 

Area{G') — ^ Area{e) ■ Weight{e). 

eeE{C') 
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Definition 38. We say that G is antirational if dG contains a connected part r which contains no subsets 
of rational slope. 

Remark 20. For A which are neither Q-polygons nor antirational, the Area(C‘^) can be infinite or finite. 
Lemma 44. If dG is antirational then Area{G^) is infinite. 

Proof. This lemma is somewhat similar to Theorem 1 in [20] but has no direct relation with it. Since is 
antirational, there is a part r of d^l which contains no sides of rational slope. We consider the polygons fig. 
Let us choose C dilg such that —>■ r. Recall that Area{G^ n fie) equals to the sum of euclidean lengths 
of the sides of with the weights, given by their directions. Since r contains no sides of rational slopes, then 
for each iV > 0 the sum Si of euclidean lengths of edges in with slopes (p, q), < N tends to zero. 

This implies that Area{rif) —>■ oo as e —>■ oo, because Area{r^) > N x {re — Se). □ 

Question: what is the euclidean length of this tropical curves? 

Since the tropical symplectic area of G^ in infinite, it is natural to consider some kind of regularisation. 

Definition 39. Regularized area of Cp is the limit by compacts iF C of Area{Gp n K) — Area{C^ n K). 

Lemma 45. If did contains no sides with rational slopes then the regularised area of each Cp is zero. 

Proof. Indeed, if we consider the weights on dfle for Cp and C^, they differ only on finite number of intervals. 
Since 917 contains no sides with rational slopes, all these intervals tends to zero. □ 

17. Discussion 

17.1. Minimal curves and genus. Sandpile model provides us a solution of the following problem: given 
a set of points in a lattice polygon A we want to find a tropical curve of minimal symplectic tropical area, 
passing through these points. More precisely, we a looking for tropical curve, which intersect A only at 
vertices of A and the intersection of the curve and A has the minimal tropical area. 

We should mention that there is a whole family of such curves - shrinking or inflating a component of the 
complement to the curve does not change its tropical symplectic area. 

So, the question remains: what is so special about our curve among all the minimal curve. Theoremj^says 
that the sand solution minimizes the integral of the toppling function, i.e. of the polynomial which defines 
the limiting curve. In fact, this means that the limiting curve has the maximal genus among all the minimal 
curves passing through these points and under above conditions. 

17.2. Generalizations to other graphs. As the reader can notice, the fact that the dimension of the lattice 
is two was not heavily used in the proofs. So, all similar results are obtained for higher-dimensional lattices, 
where we get not tropical curves but tropical hypersurfaces. This approach can be further generalized for 
any graph with a certain map rj to Z". For example, for the Cayley graph of a group G this map rj imitates 
the map to the Cayley graph of the abelianization of G. The changes in the proofs are quite straightforward. 
The only trick is how to define the boundary. We proceed by the following way. 

We start from a graph G decorated with a collection of maps Si, 0 : 2 ,..., : G —>■ Z. We require that 

all these function Xi are discrete linear harmonic functions on G. Linear means here that for each path 
Vl,V 2 ^... ,Vn of a length n we have a uniform on n estimation — a;i(ui)| < cn. Then, we define a 

boundary by the inequalities dG = {vlJ^GjXiiv) < Cj} for some chosen set of integers Cij,Cj,i G l..m,j € 
l..n. 

17.3. Sand dynamic on tropical varieties, divisors. It is easy to see that we can produce the same 
type of problems for tropical varieties, if we have a sort of grid there, given by the affine structure and maps 
Xi,... ,Xn usually given by projections on the coordinates. The convergence results can be proven in the 
same way. 

What is an interesting aspect of the possible applications is the tropical divisors. Indeed, using relaxation 
we can understand by the sand dynamic if there exist a divisor linearly equivalent to a given one, passing 
through prescribed set of points. For that we represent this divisor as a coloured set, using the patterns 
we already discovered, then we add sand to the points pi,p2,... ,Pn and relax the obtained state. If the 
relaxation terminates, it produces the divisor with the required properties. If not, that means that such a 
divisor does not exist. 
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also we understand the dynamics, the waves propagate through walls with marked points, solve this 
non-locality problem 

Also we proved that the support of each wave is a face, contractible, if you want. 
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